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SECTION  1 


INTRODUCTION 


1.1  SUMMARY 

This  report  covers  the  work  performed  by  General  Dynamics  on  the  Active  Control  of 
Space  Structures  (ACOSS  SEVEN)  program.  The  objective  of  ACOSS  SEVEN  was  to 
extend  ancTreTTne  General  Dy nimicsr modern  control  theory  approach  Model  Error 
Sensitivity  Suppression  (MESS)  and  thus  provide  a  unified  dynamics  and  control  tech¬ 
nology  base  for  large  spaceborne  structures  requiring  stringent  line  of  sight  (LOS) 
requirements.  Additional  objectives  were  analysis  and  simulation  and  demonstration 
test  planning. 

The  problem  of  achieving  comprehensive  control  of  large  spaceborne  structures  is 
essentially  one  of  mode  error.  A  control  system  must  be  designed  for  a  truncated 
model  of  an  essentially  infinite  -order  system.  The  truncated  modes  fall  into  two 
classes:  known  and  unknown.  Hie  known  modes  are  easily  handled  through  the  MESS 
algorithm  by  a  simple  alteration  of  the  optimal  performance  index.  To  cope  with  the 
problem  of  model  error  caused  by  the  truncation  of  unknown  high  frequency  residual 
modes.  General  Dynamics,  during  the  course  of  the  study,  developed  the  filter 
accommodation  algorithm.  Filter  Accommodated  Optimal  Control  allows  the  design 
of  low  bandwidth  optimal  control  and  estimation  systems  such  that  the  unknown  high 
frequency  residual  modes  are  gain  stabilized.  Filter  accommodation  allows  the 
inclusion  of  a-priori  designer  specified  low  frequency  filter  dynamics  in  the  control 
loop  such  as  Tchebycheff  or  Butterworth;  the  optimal  control  law  is  designed  specific 
to  the  chosen  filter  dynamics  such  that  the  low  pass  characteristics  of  the  compensa¬ 
tor  loop  are  not  nullified.  The  algorithm  allows  filter  attenuation  to  begin  within  the 
control  system  bandwidth,  the  controller  essentially  works  through  the  filter,  and 
low-order  filters  may  be  used  to  achieve  the  required  degree  of  residual  mode 
attenuation.  Filter  accommodation  is  compatible  with  the  basic  MESS  algorithm. 
Simultaneous  use  of  both  algorithms  -  Filter  Accommodated  Model  Error  Sensitivity 
Suppression  (FAMESS)  provides  simultaneous  alleviation  of  the  t  .-incation  problems 
of  LSS  control;  i.  e. ,  the  suppression  and  attenuation  of  model  error  caused  by  the 
truncation  of  known  and  unknown  modes. 

In  the  area  of  overall  system  stability  a  Liapunov  stability  analysis  was  performed 
for  the  filter  accommodated  designs.  A  stability  criteria  was  developed  which  is  a 
function  of  residual  mode  parameters  and  compensator  gain  levels. 

The  Draper  Evaluation  Structure  #2  was  analyzed  and  disturbance  accommodation 
techniques  in  conjunction  with  decoupled  control  were  applied  to  meet  the  LOS  error ' 
specifications. 
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As  part  of  the  Demonstration  Planning  task,  various  large  space  structure  conceptual 
designs  were  surveyed  to  Identify  the  types  of  structural  element  used.  These 
elements  were  then  put  into  a  versatile  test  structure  design  which  includes  provision 
to  vary  structural  stiffness. 

1.2  RE  COMMENDATIONS 

The  MESS/FAMESS  control  technique  has  demonstrated  significant  progress  toward 
the  goal  of  a  generic,  unified  control  methodology  for  large  spaceborne  structures 
with  stringent  stabilization  and  control  requirements:  The  truncation  and  attendant 
spillover  problems  for  LSS  control  by  Linear  Quadratic  Gaussian  (L/QG)  methods  are 
essentially  resolved.  As  a  result  of  our  ACOSS  work,  one  can  now  routinely  design 
finite  dimensional  LQG  controllers  for  infinite  dimensional  systems  employing  a 
systematic  methodology. 

However  some  additional  progress  is  required  before  the  control  theory  requirements 
can  be  declared  fulfilled  for  all  classes  of  systems  operating  in  all  disturbance 
environments.  The  resolution  of  the  fundamental  modal  truncation  problems  permit 
attention  to  be  directed  toward  measures  of  system  performance  and  merit;  i.  e. , 
gain  margins,  phase  margins,  system  sensitivity,  robustness,  disturbance  rejection 
and  model  reduction.  It  makes  little  engineering  sense  to  expend  a  large  amount  of 
time  and  effort  on  these  questions  before  the  basic  issues  of  modal  truncation  are 
formally  resolved.  In  addition  to  the  control  issues  just  enumerated,  as  technology 
moves  toward  actual  hardware  implementation  a  host  of  design  interface  and  integra¬ 
tion  problems  appear;  i.e. ,  just  how  does  one  distribute  the  damping  requirements 
among  the  structural  design  process,  passive  damping  devices,  and  the  active  control 
system  in  a  systematic  procedure  such  that  the  overall  system  is  truly  integrated  and 
no  damping  technique  is  considered  in  ad  hoc  fashion.  The  following  areas  are 
recommended  for  additional  study. 

1. 2. 1  THEORY  EXTENSION. 

Stability  Margins:  The  Liapunov  stability  results  obtained  thus  far  are  functions  of 
the  residual  mode  parameters.  One  obtains  stability  results  as  functions  of  worst 
case  parameters.  This  approach,  while  sufficient  for  some  systems,  may  be 
inapplicable  or  cumbersome  for  others.  Thus  some  attempt  should  be  made  to 
develop  stability  margins  that  are  not  dependent  on  unknown  modal  parameters. 

In  addition,  now  that  the  truncation  problems  are  resolved,  classical  gain  and  phase 
margin  concepts  should  be  generalized  to  multivariable  control  and  applied  to  the  LSS 
control  problem.  Many  such  results  exist  as  Isolated  facts  in  the  literature,  but  as 
usual,  these  results  must  be  modified  and  melded  together  to  provide  a  useful  design 
methodology. 
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System  Robustness;  The  closed-loop  robustness  properties  of  LQG  designs  Is 
ultimately  tied  to  the  optimal  compensator.  General  Dynamics  made  significant 
progress  on  the  compensator  question  via  its  constrained  compensator  algorithm 
during  ACOSS  Phase  la;  however,  recently  developed  singular  value  robustness 
measures  should  be  modified,  as  necessary,  and  Incorporated  into  the  FAMESS 
design  methodology  such  that  the  compensator  which  emerges  from  the  design 
process  has  a  prescribed  degree  of  robustness  as  measured  by  the  <ringni«i»  value 
criterion. 

Disturbance  Rejection:  Disturbance  rejection  techniques  must  be  suitably  modified 
to  be  useful  for  LSS  control.  Disturbance  site  estimation  and  feedforward  control 
is  recommended  to  cancel  disturbance  effects  whenever  possible.  This  preference 
is  for  two  reasons:  first,  the  finite  propagation  time  entailed  by  the  distributed 
nature  of  large  space  systems  could  permit  the  disturbance  process  to  be  operational 
for  long  periods  of  time  before  it  was  detected  at  nominal  system  outputs .  Estima¬ 
tion  of  the  disturbance  process  at  the  disturbance  site  eliminates  this  possibility. 
Secondly,  the  feedforward  control  technique  provides  a  two  degree  of  freedom  control 
system  such  that  the  disturbance  accommodating  loop  can  be  designed  indapandAntiy 
of  the  modal  damping  loop. 

The  Draper  Evaluation  Structure  Number  2  provides  an  outstanding  example  of  the 
need  for  the  development  of  disturbance  accommodation  technology  suitable  for  tas 
control.  No  amount  of  modal  damping  will  permit  the  system  to  meet  LOS  require¬ 
ments.  Thus  one  must  develop  additional  control  technology  which  is  compatible  with 
modal  damping  technique  of  the  active  control  system.  It  is  again  noted  that  the 
development  of  these  secondary  control  technologies  for  LSS  control  was  impossible 
until  the  truncation  problems  were  resolved. 

Closed- Loop  Model  Reduction:  Some  means  of  controller  simplicatlon  should  be 
attempted;  i.  e. ,  one  first  designs  the  high  order  controller  and  then  systematically 
eliminates  controlled  states  and  readjusts  system  gains  to  maintain  performance. 

The  end  result  of  this  process  is  a  high  performance  controller  of  reduced  complexity. 

1.2.2  ROLE  OF  DISTURBANCE  ACCOMMODATION 

The  Draper  Evaluation  Structure  #2  raises  the  question  of  the  relative  roles  of  active 
structural  damping  and  disturbance  accommodation.  Although  active  structural  damp¬ 
ing  technology  has  received  by  far  most  of  the  emphasis  in  prior  active  structural 
control  work,  the  Draper  Evaluation  Structure  #2  could  not  be  made  sufficiently  quiet 
to  meet  specifications  by  the  application  of  active  damping.  Thus  it  is  recommended 
that  studies  be  performed  to  establish  whether  disturbance  accommodation  technology 
should  receive  increased  emphasis  or  whether  the  Evaluation  Structure  presents  a 
problem  which  is  best  solved  by  structural  design  changes. 
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1.2.3  STRUCTURAL  CONTROL  COMPONENTS 

Reviews  of  the  active  structural  control  state  of  the  art  clearly  indicate  that  control 
theory  is  far  ahead  of  practical  implementation.  Except  for  a  few  cases,  active 
structural  control  test  work  has  used  actuators  and/or  sensors  which  are  not  appli¬ 
cable  to  use  on  real  structures  in  space.  Thus  additional  theory  developments  are 
certainly  important  but  increased  emphasis  on  practical  control  system  implementa¬ 
tion  appears  to  be  essential  if  the  active  structural  control  technology  is  to  be  ready 
to  support  large  space  systems  in  a  timely  manner. 
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SECTION  2 


CONTROL  THEORY 


A  systematic  procedure  for  handling  the  truncation  of  unknown  high  frequency  re¬ 
sidual  modes  has  long  been  an  open  problem  in  Large  Space  System  (LSS)  control, 
(References  1-6).  Residual  modes  are  those  vibration  modes  that  are  neglected  during 
the  design  process.  These  modes  are  usually  poorly  known  and  of  high  frequency. 

(The  precision  of  modal  calculations  decreases  as  modal  frequency  increases. )  This 
problem  of  unknown  dynamic  characteristics  at  high  frequencies  is  endemic  throughout 
control  engineering  and  is  not  generic  to  the  LSS  control  problem  specifically.  The 
designer  rarely  knows  the  high  frequency  characteristics  of  the  plant;  however,  due  to 
the  infinite  series,  mods’  representation  of  large  space  systems,  the  problem  of  un¬ 
known  high  frequency  m<  .ai  dynamics  is  more  visible  and  plays  a  greater  part  in  the 
LSS  design  process  than  in  many  other  control  systems. 

Classically,  the  usual  technique  to  cope  with  unknown  dynamics  at  high  frequencies 
has  been  the  roll-off  filter,  i.e. ,  loop  gain  decreases  as  frequency  increases,  and 
thus  all  dynamics  above  a  certain  frequency  are  gain  stabilized,  (Reference  7).  This 
is  the  approach  that  we  shall  adopt  here;  however,  this  roll-off  filter  concept  is  not 
easily  applied  in  optimal  control  and  estimation  systems,  which  are  generally  high 
bandwidth  systems.  Thus,  we  develop  a  new  specific  optimal  control  algorithm  that 
permits  the  inclusion  of  designer  selected  filter  dynamics.  The  primary  reason  for 
this  roll-off  filter  approach  stems  from  the  uncertainty  of  LSS  plant  dynamics  at  high 
frequencies.  Ultimately,  we  know  nothing  about  the  high  frequency  modes;  neither 
mode  shape  nor  frequency.  Such  lack  of  knowledge  precludes  a  precise  mathematical 
formulation  of  the  problem.  (Recall  that  characteristics  of  sensors  and  actuators, 
both  physical  devices,  are  not  known  with  great  precision  at  high  frequencies. ) 

Thus  we  adopt  the  roll-off  filter  approach  and  attenuate  what  we  do  not  know.  How¬ 
ever,  this  roll-off  filter  approach  must  be  further  modified  for  LSS  control  because 
of  the  extremely  dense  modal  spectrum  of  the  LSS.  As  illustrated  in  Figure  2-1,  the 
previous  approach  for  small  satellites  was  the  design  of  low  bandwidth  controllers 
and  a  stiff  structure  such  that  all  flexible  body  modes  were  outside  the  control  system 
bandwidth,  but  this  approach  is  not  feasible  for  large  systems  where  the  vibration 
modes  are  fractions  of  a  Hertz  and  frequency  separations  are  on  the  order  of  a  one- 
one  hundredth  of  a  Hertz.  Thus  we  must  not  only  attenuate  residual  modal  response, 
but  must  actually  design  a  controller  that  works  through  the  filtering  system  such 
that  certain  modes  deemed  critical  to  the  control  task  are  controlled  and  the  responses 
from  unknown  residual  modes  are  attenuated. 
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Figure  2-1.  Classical  Cldsed-Loop  Control  Approach 
2. 1  THE  LSS  CONTROL  PROBLEM 

The  LSS  control  problem  can  be  sufficiently  described  by  recourse  to  three  state 
space  spacecraft  models:  the  evaluation  model,  the  design  model,  and  the  control 
model,  (Reference  8). 

2. 1. 1  EVALUATION  MODEL.  This  is  the  most  complex  model  and  has  the  closest 
fidelity  to  the  actual  physical  system.  This  model  is  used  to  validate  the  closed-loop 
system  after  design  and  usually  encompasses  unknown  modes,  disturbances,  and  non- 
linearities.  In  actuality,  the  evaluation  model  can  be  any  set  of  dynamic  equations 
that  can  be  placed  on  a  computer.  Figure  2-2  characterizes  a  typical  linear  high- 
order  evaluation  model  that  is  employed  to  evaluate  truncation  error  sensitivity. 


The  following  definitions  apply  to  the  pardoned  equations: 


Controlled  states:  x.  eRn 

C 

Suppressed  states:  xs  «RS 
Residual  states:  xr  «Rr 
Control  vector:  u  e  Rm 
Observation  vector:  y  eR* 

Control  spillover:  Beu  and  Bru 
Observation  spillover:  Csxg  and  C^Xj, 


Xc 

Ac 

0 

0 

xc 

Be 

*s 

= 

0 

As 

0 

Xs 

+ 

BS 

.xr. 

0 

0 

Ar. 

LBrj 

[y]  ~  [Cc  Cs  Crj 


xc 

xs 

xr 


The  residual  state  vector  xr  is  of  arbitrary  dimension. 


Figure  2-2.  Evaluation  Model 


The  plant  matrix  has  the  block  diagonal  form  shown  because  the  state  variables 
characterize  modal  coordinates. 

This  partition  of  the  state  vector  reflects  control  objectives  and  modeling  accuracy. 

The  controlled  states,  Xc,  characterize  those  critical  modes  that  must  be  controlled 
to  achieve  satisfactory  system  performance.  The  suppressed  states,  xs,  characterize 
those  modes  whose  control  is  not  critical  to  the  control  task.  The  suppressed  modes 
are  well  known,  but  do  not  greatly  influence  the  control  objective.  It  is  desirable  to 
eliminate  these  modes  from  the  dynamic  model  in  order  to  minimize  system  complex¬ 
ity.  The  residual  states,  xr,  characterize  unknown  or  unreliably  computed  modes. 

The  residual  modes  are  usually  of  high  frequency  and  are  neglected  during  the  design 
process.  These  residual  modes,  however  poorly  known,  do  constitute  system  charac¬ 
teristics  with  which  the  design  process  must  cope.  It  is  these  residual  states  that  will 
receive  primary  emphasis  in  this  development. 

The  spillover  terms  indicate  undesirable  effects  of  reduced-order  control  system 
design.  Control  spillover  is  the  excitation  of  uncontrolled  modes  by  the  controller, 
and  observation  spillover  is  the  sensing  of  uncontrolled  modal  responses  by  the 
observer.  In  a  closed-loop  system,  simultaneous  control  and  observation  spillover 
can  cause  instability,  (References  5,  6).  Thus,  spillover  reduction  is  one  objective 
for  low-order  control  of  high-order  systems. 

2. 1. 2  DESIGN  MODEL.  This  model  contains  states  corresponding  to  those  modes 
that  can  be  reliably  computed  and  can  be  viewed  as  a  truncated  version  of  the  evalua¬ 
tion  model.  Figure  2-2,  with  the  residual  states  deleted.  Figure  2-3  characterizes 
the  typical  linear  design  model.  The  design  model  usually  contains  too  many  states  to 
be  actively  controlled.  Also,  many  of  the  states  labeled  xg  are  not  critical  to  the 
performance  objectives;  thus,  the  design  model  is  further  truncated. 

2.1.3  CONTROL  MODEL.  This  model  contains  those  states  that  correspond  to  critical 
modes  that  must  be  controlled  to  meet  performance  requirements.  Figure  2-4 
characterizes  the  typical  linear  control  model. 

If  an  identity  estimator  is  employed  in  the  system  design,  the  low-order  control  model 
yields  a  model  for  the  estimator. 

2.1.4  LSS  CONTROL  PROBLEM.  The  LSS  control  problem  now  can  be  cast  simply 
as:  design  a  controller  that  stabilizes  the  control  model.  Figure  2-4,  remains  stable, 
and  meets  performance  criteria  when  applied  to  the  evaluation  model.  Figure  2-2. 

Alternately  we  pose  the  question:  is  it  possible  to  design  a  controller  based  on  the 
dynamic  system  Xg,  and  suppress  sensitivity  to  modeling  errors  deliberately 
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where  xc  is  the  state  vector  to  be  controlled 

xs  is  the  state  vector  to  be  suppressed 
y  is  the  measurement  vector  based  on  xc  &  xs 
Ac.  Be,  Cc,  As,  Bs,  Cs  are  constant  matrices 
vc.  vs,  w  are  white  Gaussian  noise 
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Figure  2-3.  Design  Model 


Xc  =  AcXc  +  Bcu  +  v 
yc  =  CqXc  +  w 


where  xc  is  state  vector  to  be  controlled 
u  is  control  vector 
yc  is  measurement  based  on  xc 
Ac,  Bc,  Cc  are  constant  matrices 
v,  w  are  white,  Gaussian  noise 
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Figure  2-4.  Control  Model 


introduced  by  truncating  states  Xg  and  xr.  Thus  our  truncation  problem  is  twofold,  we 
must  deal  with  the  truncation  of  known  and  unknown  modal  states. 

The  basic  Model  Error  Sensitivity  Suppression  (MESS)  algorithm  copes  with  the 
truncation  of  known  modal  characteristics  by  incorporating  information  from  the  de¬ 
sign  model  in  the  optimal  performance  index  such  that  the  states  xs  receive  reduced 
controller  excitation,  (References  8-13).  These  known  suppressed  states  are 
essentially  decoupled  from  controller  action. 

However,  complete  solution  to  the  LSS  control  problem  requires  some  means  to  cope 
with  the  truncation  of  the  unknown  residual  states  xr  of  the  evaluation  model.  Recall 
that  we  have  essentially  zero  knowledge  concerning  these  modal  states,  except  perhaps 
location  of  a  lower  bound  in  the  modal  frequency  spectrum.  Thus  we  adopt  the  roll-off 
filter  approach  and  attenuate  what  we  do  not  know,  but  with  the  added  stipulation  that 
the  control  algorithm  must  be  capable  of  working  through  the  filter. 

The  outline  of  this  report  is  as  follows: 

1)  We  develop  a  filter  accommodation  algorithm  that  is  compatible  with  optimal 
control  and  Model  Error  Sensitivity  Suppression  (MESS).  Filter  accommodation 
allows  the  designer  to  select  the  roll-off  filter  dynamics  as  required  and  de¬ 
velops  an  optimal  control  algorithm  specific  to  the  chosen  filter  dynamics.  Thus 
Filter  Accommodated  Model  Error  Sensitivity  Suppression  (FAMESS)  provides 
complete  solution  for  the  truncation  of  known  and  unknown  modal  coordinates. 

2)  A  simple  two -mode  design  model  using  second  order  filters  is  presented. 

3)  Extensions  to  the  filter  accommodation  technique  are  briefly  outlined.  These 
include  output  feedback,  and  through  duality  theory  Kalman  filtering.  Our  re¬ 
sults  are  preliminary  on  this  last  item. 

2. 2  FREQUENCY  RESPONSE  AND  FULL  STATE  FEEDBACK 

One  of  the  characteristics  of  closed-loop  systems  designed  by  optimal  control  methods 
is  high  bandwidth.  In  many  systems  the  optimality  condition  tends  to  undo  (cancel 
might  be  more  descriptive)  all  dynamics  and  rearrange  parameters  such  that  the  re¬ 
sulting  system  has  a  flat  frequency  response  in  the  controlled  frequency  region  and 
rolls  off  as  (1/s)  at  high  frequencies.  This  tendency  of  optimal  systems  to  be  of  high 
bandwidth  is  well  noted  in  the  literature,  (References  14-19).  This  high  bandwidth 
behavior  is  appropriate  when  system  dynamics  are  well  known,  or  when  the  physical 
system  is  relatively  noise-free.  Unfortunately,  such  conditions  are  rarely  met  in 


practice.  Indeed,  one  of  the  reasons  that  optimal  regulator  techniques  have  gained 
so  little  acceptance  in  the  industry  (Kalman  filters  excepted)  has  been  the  inability  to 
translate  system  bandwidth  requirements  into  adequate  performance  indices. 

An  analysis  of  the  roll-off  problem  reveals  that  the  full  state  feedback  system,  optimal 
or  non-optimal,  results  in  an  open-loop  gain  matrix,  L0(s),  (the  product  of  the  plant 
and  compensator  matrices)  which  usually  rolls  off  as  (1/s).  Thus,  we  see  that  as  the 
optimal  regulator  belongs  to  the  class  of  full  state  feedback  systems, -it  usually  be¬ 
haves  as  (1/s)  at  high  frequencies. 

Thus,  the  solution  to  the  roll-off  problem  becomes  obvious:  Remove  the  optimal 
regulator  from  the  class  of  full  state  feedback  systems.  Hence,  we  develop  a  specific 
optimal  control  algorithm  that  does  not  require  complete  state  feedback  and  is  also 
compatible  with  the  MESS  algorithm. 

It  should  be  re-emphasized  that  the  MESS  algorithm  provides  solution  for  the  problem 
of  known  modal  truncation.  The  problem  we  address  here  is  not  directly  amenable  to 
these  techniques  as  the  high  frequency  dynamics  are  unknown. 

Some  pertinent  characteristics  of  full  state  feedback  systems  are  illustrated  by  the 
simple  system  of  Figure  2-5.  The  system  is  of  order  two  and  we  feed  back  two  state 
variables,  say  position  and  velocity.  Thus,  the  return  difference  which  equals  one 
plus  the  loop  gain,  L0(s),  will  have  the  structure  shown.  In  general,  as  in  the  second 
order  case,  the  numerator  order  of  the  entries  LQ(s)  is  one  less  than  that  of  the  de¬ 
nominator;  these  entries  behave  as  (l/s)  at  high  frequencies.  It  is  this  property  of  full 
state  feedback  systems  that  we  seek  to  remove  from  the  optimal  regulator.  Thus  our 
basic  approach  is  to  find  an  optimal  system  that  has  gains  for  certain  pre-selected 
filter  states  which  are  equal  to  zero.  Essentially,  we  are  seeking  an  optimal  solution 
to  the  regulator  problem  that  does  not  undo  our  filter  dynamics.  The  primary 
property  that  we  wish  to  preserve  at  this  point  is  the  high  frequency  attenuation 
characteristics  of  our  filters.  Hence,  we  resort  to  specific  optimal  control  in  order 
to  preserve  the  desirable  properties  of  optimal  design,  and  simultaneously  retain  our 
control  of  the  system  bandwidth. 

2.3  CONCEPTUAL  BASIS  FOR  THE  FILTER  ACCOMMODATION  ALGORITHM 

Central  to  the  filter  accommodation  algorithm  is  limited  state  feedback  from  a  tan- 
demly  connected  system.  As  we  have  seen,  full  state  feedback  results  in  first-order 
high  frequency  behavior.  In  order  to  achieve  the  roll-off  desired  we  generate  the 
system  structure  shown  in  Figure  2-6.  The  state  feedback  gains  around  the  plant 
(states  xj)  generally  result  in  a  first-order  roll  off.  The  filter  is  assumed  to  be  of 
order  f.  This  system  transfer  function  matrix  from  r  to  W*U  roll,  off  at  least  as 
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Figure  2-5.  Return-Difference  for  a  Second-Order  System 

fast  as  (1/s)  *  +  Thus,  the  designer  can  specify  the  high  frequency  system  be¬ 
havior  desired. 

Now  we  desire  the  gain  matrix  K  =  (Kj  O]  to  be  optimal.  This  poses  a  specific 
optimal  control  problem:  minimize  an  appropriate  performance  index,  simultaneously 
stabilize  system  I,  and  generate  zero  gains  for  system  II. 

Specific  Optimal  Control  (S.O.C.)  was  developed  circa  1965,  pre-Luenberger 
Observer,  to  cope  with  the  output  feedback  problem.  The  original  technique  ("Optimal 
Control  with  Unavailable  States,"  John  F.  Cassidy,  Ph.D. ,  Thesis,  R.P. I. ,  1969: 
Reference  20)  was  iterative  in  nature  and  required  the  solution  of  a  sequence  of 
optimal  regulator  problems.  Sequential  adjustments  to  the  weigting  matrices  achieved 
the  desired  results. 

The  technique  we  develop  here  differs  in  scope  and  solution  procedure:  1)  It  is  not 
iterative,  and,  2)  it  is  not  an  output  feedback  scheme,  but  rather  an  "internal"  feed¬ 
back  problem  where  state's,  xj  behind  the  filter  subsystem  are  fed  back.  Thus,  in  the 
final  implementation,  an  estimator  will  be  employed,  but  in  our  control  law  we  will 
only  feed  back  estimates  of  the  substate  vector  x^.  If  the  filter  states,  Xg,  are  easily 
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Figure  2-6.  Conceptual  Diagram  for  Filter  Accommodation 

accessible,  reduced  order  estimation  can  be  employed,  i.e. ,  in  many  situations  we 
would  not  estimate  the  filter  states.  Hence,  our  basic  approach  is  to  develop  state 
gain  matrices  of  the  form  shown  in  Figure  2-6  and  feed  back  estimates  of  these 
states  as  required. 

As  Figure  2-7  illustrates,  the  closed-loop  system  structure  that  results  from  appli¬ 
cation  of  the  filter  accommodation  technique  consists  of  the  tandem  connection  of 
plant,  filters,  ami  optimal  compensator.  An  additional  feedback  loop  from  compen¬ 
sator  to  the  filtering  system  results  from  application  of  the  algorithm.  This  second¬ 
ary  feedback  path  is  static  in  nature  (a  matrix  of  constant  gains)  and  compensates  for 
the  effect  of  zero  feedback  gains  from  the  filter  states.  The  filter  accommodation 
matrix  also  functions  as  an  additional  input  matrix  for  the  filtering  system. 

2.4  OVERVIEW  OF  THE  ALGORITHM  DERIVATION 

The  filter  accommodation  algorithm  is  developed  as  follows: 


2-9 


(ESTIMATOR-BASED) 


150S1275-13 


Figure  2-7.  Control  System  Structure  for  Prespecified  Filter  Dynamics 

1)  Formulate  the  problem  as  that  of  a  -tandemly  connected  system  such  that  a  first 
dynamic  system  drives  a  second  dynamic  system,  but  not  vice  versa. 

2)  Form  the  algebraic  Riccati  equation  for  the  tandem  system.  This  allows  us  to 
observe  the  optimal  coupling  which  would  exist  in  a  full  state  feedback  optimal 
control  solution. 

3)  Impose  the  zero  gain  conditions  on  the  partitioned  Riccati  equation.  This  allows 
us  to  look  for  simplification  and  possible  solution  techniques. 

4)  Observe  that  after  imposition  of  the  zero  gain  conditions,  the  algebraic  Riccati 
equation  is  decoupled  into  three  subsystems  that  allow  sequential  solution:  A 
Liapunov  equation  for  the  filter  (system  II);  a  Sylvester  equation  that  couples  the 
plant  (system  I)  with  the  filters  (system  II).  And  a  modified  Riccati  equation 
which  provides  stabilization  for  system  I  (the  plant). 
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Figure  2-8  provides  the  detailed  view  of  the  tandem  system  which  illustrates  the  nota¬ 
tion  and  symbology  that  we  shall  use  in  our  development  of  the  filter  accommodation 
algo:,vhm. 
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Figure  2-8.  Tandem  System  Structure 

System  (A.^,  corresponds  to  the  plant  and  System  (Agg)  corresponds  to  the  filters. 
Matrix  A2y  corresponds  to  the  coupling  between  the  output  of  the  plait  amt  the  input  to” 
the  filtering  system.  The  matrix,  B2p  indicated  by  the  dashed  lines,  initially  zero, 
is  determined  during  the  design  process.  The  matrix  B21  functions  as  an  additional 
input  matrix  for  the  filter  system.  B2j  is  determined  so  as  to  force  the  zero-gain 
solution  in  the  state  regulator.  The  introduction  of  matrix  B2j  eliminates  the  need 
for  iteration  present  in  Cassidy's  method. 

Matrix  representation  of  the  tandem  system  yields  the  form  shown  in  Figure  2-8.  The 
zero  submatrix  in  the  upper  right-hand  corner  of  the  A -matrix  implies  the  tandem 
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connection:  Matrix  A12  is  identically  equal  to  zero;  thus,  system  states  xj  drive 
system  states  x2,  but  states  x2  do  not  in  turn  drive  states  xj.  The  coupling  is  only 
unidirectional,  as  shown  in  the  block  diagram. 

2.  5  DERIVATION  OF  THE  ALGORITHM 

Direct  substitution  of  the  tandem  equations  into  the  algebraic  Riccati  equation  yields 
the  coupled  algebraic  equations  shown  in  Figure  2-9.  Our  primary  emphasis  will  be 
on  the  quadratic  term  involving  R -inverse,  and  we  will  decouple  these  equations  with 
an  appropriate  choice  of  the  B2j  matrix. 

To  simplify  unwieldy  algebraic  expressions,  we  define'the  matrices  W^,  as 

shown  in  Figure  2-10,  and  indicate  the  term  by  term  expansion  of  the  quadratic 
factors.  In  sequel  we  shall  set  W,,^  identically  equal  to  zero  and  thus  decouple  our 
equations. 

With  reference  to  Figure  2-11,  analysis  of  the  optimal  feedback  gain  matrix  provides 
the  constraint  condition  that  we  desire:  zero  feedback  for  states  x2«  Thus,  following 
equation  (1)  we  partition  the  gain  matrix  into  submatrices  Kj  and  Kg.  We  next  apply 
the  relationship  for  optimal  feedback  control  (equations  2  and  3).  Expression  of  the 
optimal  control  in  terms  of  the  W-matrices  yields  equation  (4).  Imposition  of  our 
constraint  condition  yields  equation  (5).  Thus,  we  merely  solve  equation  (5)  for  B21  , 
as  shown  in  equation  (6).  In  sequel  we  show  that  P22- inverse  always  exists,  thus 
permitting  solution  for  B21. 

We  note  that  our  solution  is  not  the  only  one  which  drives  K ^  to  zero.  If  in  equation 
(3),  B2i  is  set  equal  to  zero,  the  constraint  relationship  becomes  Bj2T  P12  =  O. 

Thus,  an  alternate  solution  would  require  P^2  to  lie  in  the  kernel  of  B^T.  At  this 
time,  it  is  thought  that  such  a  solution  may  be  too  restrictive;  however,  this  null 
space  solution  is  an  interesting  facet  of  the  algorithm  in  that  it  would  eliminate  the 
additional  feedback  loop  required  for  B2j.  Such  a  development  may  be  attractive  for 
certain  classes  of  control  problems. 

The  previous  development  provided  the  motivation  and  conceptual  overview  for  the 
filter  accommodation  technique.  We  analyzed  the  roll-off  problem  and  decided  full- 
state  feedback  was  not  desirable.  We  then  translated  the  roll-off  problem  to  a  specific 
optimal  control  problem  such  that  gains  from  the  filter  states  were  zero.  Following 
this,  the  conditions  for  zero  filter  state  gains  were  developed.  It  yet  remains  to  be 
demonstrated  that  our  constraint  set  yields  non-conflicting  solvable  equations. 


2-12 


<*11*11  +  *12*21) 
<*12*11  +  *22*21 


l"'1!  H*11*11  +  *{l  >12)  <*11  *12  +  *Jl*22» 


«11 

*21 


I"'1)  [wj. 


M 


*11 

iR-niwliwiii 

MllR'Iwl,  Vul'Ml* 

*2i 

W*,R-’wli  W2,R-1wli 

TERM-BY-TERM  EXPANSION  VELDS 

WllR*%?l  =  1*11*11  +  *12*2ll  l*’1)  1*11*!  1  +  *21*12) 

WziR-’W^  =  1*^*1,  +  *22*21)  l**1!  1*11*11  +  *21*121 

*11*' **21  =  1*11*11  +  *12*21)  l*'1)  1*11*12  +  *21*22) 

WjlR-1*!,  *  1*12*11  +  *22*2lJ  I*1)  1*11*12  +  *21*22J 

Figure  2-10.  Expansion  of  the  Quadratic  Riccati  Term 
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Figure  2-11.  Analysis  of  the  Optimal  Feedback  Gain  Matrix 
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Figure  2-12.  Zero-Gain  Condition  and  the  Quadratic  Riccati  Term 

As  noted  in  Figure  2-12,  the  gain  matrix  is  partitioned  into  two  submatrices,  Kj  and 
Kg,  and  can  be  expressed  in  terms  of  and  W21,  as  shown  in  equation  (1).  The 
condition  Kg  =  O  implies  Wgj  =  O,  and  all  quadratic  terms  involving  W21  are  zero 
matrices,  as  shown  in  equations  (3-5).  Thus  the  quadratic  term  assumes  the  form 
shown  in  equation  (6),  and  the  quadratic  submatrix  involving  W-q  can  be  expressed  as 
in  equation  (7).  The  sparse  nature  of  the  quadratic  matrix,  equation  (6),  ultimately 
results  in  decoupled  equations  allowing  sequential  solution. 

Figure  2-13  continues  the  development*.  The  partitioned  Riccati  equation  for  the  tan¬ 
dem  system  with  zero  gain  constraints  imposed,  takes  the  form  shown  in  equation  (1). 
Equation  (1)  is  concise  notation  for  four  matrix  equations  which  we  indicate  symboli¬ 
cally  as  equation  (11),  equation  (12),  equation  (21),  and  equation  (22).  (The  notation 
corresponds  to  that  of  matrix  analysis.)  In  actuality,  we  need  only  analyze  three 
matrix  equations  as  (12)  and  (21)  are  transposes  of  each  other  because  P  is  a  sym¬ 
metric  matrix. 

Expansion  of  the  four  matrix  equations  yields  the  relations  shown  in  Figure  2-14. 
Equation  (11)  is  a  modified  Riccati  equation.  Equations  (12)  and  (21)  are  Sylvester's 
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+  *21p12>  (*11p12  +  *21p22)  +  p12*2l)  p12*22 
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Figure  2-13.  The  Riccatl  Equation  with  Zero-Gain  Conditions  Imposed 
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EQUATION  (11)  -  MODIFIED  RICCATI  EQUATION 

(AVu  +  AT2iPi2)  +  (PllAn  +  Pi2*2l)  *  <Nl»11«"MlPl1  + 


+  P1282iN-1B?1Pii  -  Pl2B2lR'1B21Pi2)  +  Oil  =  0 


EQUATION  (12)  -  SYLVESTERS  EQUATION 

Al1p12  +  Pi  2*22  +  *21p22  +  0l2  =  0  (2) 

EQUATION  (21)  -  SYLVESTER’S  EQUATION  [EON’S  (12)  AND  (21)  ARE  TRANSPOSES] 
p12Al1  +  A22P12  +  P22A2T >Qi2  =  0‘T  <3) 

EQUATION  (22)  -  LIAPUNOV  EQUATION 

*22p22  +  P22*22  +  O22  =  0  ^ 
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Figure  2-14.  Zero-Gain  Riccati  Subsystems 

equations.  Equation  (22)  is  a  Liapunov  equation.  Thus  all  matrix  equations  are  well 
known,  and  are  solvable  by  currently  available  computer  algorithms  (Reference  21). 

We  next  analyze  the  three  matrix  equations  from  the  viewpoint  of  sequential  solvability 
in  an  overall  control  algorithm.  The  solution  procedure  is  outlined  in  Figure  2-15. 

We  first  analyze  the  Liapunov  equation. 

The  Liapunov  equation  is  decoupled  from  the  remaining  equations.  Recall  that  matrix 
A22  corresponds  to  the  filter  dynamics.  As  the  filtering  system  is  assumed  to  be 
stable  and  under  the  designer's  control,  this  equation  always  has  a  solution.  Thus 
P22  Is  a  positive  definite  symmetric  matrix,  and  P22-inverse  always  exists.  The 
weighting  matrix,  Q22  Is  free  and  can  be  used  to  condition  the  solution  for  P22, 
i.e. ,  increase  or  decrease  its  norm.  However,  the  fact  that  Liapunov’s  equation 
is  not  a  function  of  the  remaining  submatrices  is  Important  to  our  development  be¬ 
cause  it  permits  sequential  solution  of  the  optimal  control  problem. 

Analysis  of  Sylvester's  equation  reveals  that  it  is  a  function  of  P22,  but  not  of  any 
remaining,  and  as  yet  unknown,  submatrices.  Matrix  P22  is  known,  and  matrixQ12 
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1)  SOLVE  THE  LIAPUNOV  EM/ATVN  FOR  P22 
*22*22 +  *22*2* +  ***  =  • 

*)  SOLVE  SYLVESTER'S  EQUATION  FOR  P12 

*11*12  +  *12*22  =  “(*12  +  *21*22)  <2> 

3)  SOLVE  THE  FILTER  ACCOMMODATION  EQUATION  FOR  S21 

*21  =  -*22*12*11  <3> 

4)  SOLVE  THE  MODIFIED  RICCATI  EQUATION  FOR  Pu 

(*11-Mh)t*ii  +  PnlAn-Mn)  +  P^B^R-lR^  +  0,,  =  0  (4) 

WHERE 

■*11  =  *1l"-1*21*12  <«*> 

*11  =*11  +(*21*12+*12*2l)-*12*21*‘1*21*11  2*  (4b) 

4*)  ALPHA-SMFT  TECHNIOUES  CAN  BE  USED  TO  MODIFY  Mn 
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Figure  2-15.  Algorithm  Solution  Steps 

is  arbitrary;  thus,  we  can  solve  Sylvester's  equation  for  This  equation  has  solu¬ 

tion,  provided  A-q  and  A22  have  no  common  eigenvalues.  As  Aj^  is  generated  by  the 
plant  and  A22  is  generated  by  the  filters,  this  condition  is  easily  met  and  imposes  no 
severe  restriction  on  the  solution  process.  At  this  stage,  enough  information  exists 
to  permit  calculation  of  B2j,  the  filter  accommodation  matrix,  which  satisfies  B21  = 
-p22  P22T  ®11*  We  also  note  that  the  designer  has  a  large  degree  of  freedom  in 
shaping  B21,  because  P22  and  P12  are  nearly  free.  An  alternate  solution  procedure  is 
to  choose  Pg2  and  Pj2  such  that  they  satisfy  their  respective  equations  and  treat  Q22 
and  Q12  as  slack  variables,  which  are  determined  by  choice  of  P22  and 

We  can  now  solve  the  modified  Riccati  equation.  All  submatrices  are  known  except 
Pjl  which  is  the  solution  of  the  modified  Riccati  equation.  The  Riccati  equation  may 
be  simplified,  as  shown  in  equation  (4),  Figure  2-15. 

The  additional  matrix  which  augments  A-q,  functions  as  an  a  -shift  matrix  and 
hence  poses  no  problems  (Reference  15).  (M^  can  also  be  treated  as  a  parameter 
arising  from  a  generalized  performance  index  which  involves  a  cross  product  term  be¬ 
tween  the  states  and  control  vector,  but  we  don’t  adopt  that  approach  here.) 
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The  modified  state-weighting  matrix  Q-q  poses  no  problem  as  the  terms  comprising  it 
are  known  with  the  exception  of  which  is  an  arbitrary  positive  semi-definite 
matrix.  Thus  can  be  adjusted  to  ensure  that  is  a  positive  definite  (or  semi- 
definite  matrix  as  required. 

Hence,  all  equations  resulting  from  the  filter  accommodation  technique  are  solvable 
for  the  appropriate  system  gains,  as  shown  in  Figure  2-15. 

With  regard  to  Figure  2-15,  the  solution  steps  (1)  through  (4)  have  been  discussed  in 
detail.  The  additional  step  (4*)  simply  points  out  that  if  the  M^  —matrix  arising 
from  the  algorithm  results  in  an  unsatisfactory  eigenvalue  spectrum,  a  -shift 
techniques  can  be  employed  to  readjust  the  closed-loop  system  eigenvalues. 


OUR  APPROACH  IS  EQUIVALENT  TO  SOLVMG 
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Figure  2-16.  Complete  Optimal  Control  Problem 

It  should  be  noted  that  the  filter  accommodation  technique  is  compatible  with  the  MESS 
algorithm:  i.  e. ,  the  MESS  algorithm  works  through  the  performance  index  via  the 
control  weighting  matrix  R.  The  filter  accommodation  technique  does  not  adjust  R, 
but  rather  permits  zero  gain  state  feedback  solutions  via  the  filter  accommodation 
matrix  B21. 
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The  sequential  solution  of  the  opdnial  control  problem  just  presented  is  equivalent  to 
solving  the  standard  optimal  control  problem  illustrated  in  Figure  2-16.  Substitution 
of  the  state  weighting-submatrices  used  in  the  sequential  solution,  and  the  calculated 
value  for  the  initially  unknown  input-matrix  Bg^  will  result  in  the  specific  zero  gain 
solution  shown  in  equation  (3).  However,  it  must  be  emphasized  that  Bgj  is  unknown 
initially.  (It  can  be  shown  that  direct  solution  for  Bg^  involves  a  third-order  matrix 
equation. )  The  fact  that  our  sequential  solution  is  optimal  is  important  because  all 
the  qualities  of  optimality,  which  do  not  derive  from  full  state  feedback,  apply  to  our 
closed-loop  system.  Thus,  although  during  the  derivation  and  solution  process  we 
appear  to  have  strayed  from  "optimality"  this  is  not  the  case  and  we  do  indeed  have  a 
solution  to  the  optimal  regulator  problem  that  requires  zero  feedback  of  the  filter 
states.  This  was  our  original  objective.  The  additional  input  matrix  Bg^  for  the  filter 
states  that  is  constructed  during  the  process  provides  an  additional  degree  of  freedom 
such  that  we  can  obtain  the  required  zero  gains  in  our  feedback  matrix. 

The  final  closed-loop  system  with  estimator  in  place  takes  the  form  shown  in  Figure 
2-17  where  (Ap,  Bp,  C  )  correspond  to  the  plant  (system  I  of  our  derivation)  and 
(Af,  Bf,  Cf)  correspond  to  the  filtering  system  (system  II  of  our  derivation).  The 
matrix  product  BfCp  corresponds  to  the  coupling  matrix  Ag^:  (with  respect  to  our  de¬ 
rivation:  Ap  =  Ajj,  Af  =  Agg,  Cp  =  Cfj,  Bp  =  Bjj,  Bm  =  B21,  Kp  =  Kj,  xp  =  Xj, 
Xf  =  x2). 

The  estimator  is  an  identity  optimal  estimator  designed  via  duality  theory.  Note  that 
only  an  estimate  of  the  controlled  plant  states  is  fed  back  to  the  plant.  Also,  the 
matrix  Bm  provides  a  minor  feedback  loop,  from  the  estimator  output  to  the  internal 
driver  of  the  filter. 

The  filter/esttmator  system  including  the  minor  feedback  loop  can  be  thought  of  as  a 
compensator  attached  between  the  plant  output,  yp,  located  at  point  "a"  and  the  plant 
input,  u,  located  at  point  ”b".  As  usual,  our  definition  of  the  compensator  is  that 
system  that  has  the  sensor  measurements  as  its  input  and  the  control  vector  as  its 
output,  i.e. ,  the  transfer  function  matrix  between  points  "a"  and  "b"  of  Figure  2-17. 

In  order  to  distinguish  between  this  filter  accommodated  compensator,  and  the  usual 
optimal  compensator  without  the  additional  feedback  loop  through  Bm,  we  adopt  the 
following  nomenclature:  The  term  "compensator"  will  imply  the  usual  compensator  of 
textbook  optimal  control,  and  the  term  "filter-compensator"  will  imply  a  filter 
accommodated  compensator  with  a  filter  and  minor  feedback  loop  to  the  filter. 

After  analysis  of  the  control  system  structure,  we  may  form  the  matrix  representa¬ 
tion  of  the  plant-filter  combination  shown  in  Figure  2-18.  The  correspondence 
between  the  derivau'-m  matrices  and  the  plant-filter  matrices  is  indicated. 


Figure  2-17 .  Detailed  View  of  Filter  Accommodated  System 


We  held  the  previous  notation  as  an  analytical  aid;  however,  during  the  actual  design 
process  it  is  more  convenient  to  work  with  the  new  nomenclature. 


THE  STATE  EQUATIONS  FOR  THE  PLANT-FILTER  COMBMATNW 

VI  [*p  °  M  \h 

=  +  M 

if  BfCp  Af  Xf  Bm 

Iyil  =  lo  cf j  xp 

THUS 

Ap  =  An  Bp  =  B11  Cf  =  C12 

*1  =  *22  8m  =  B21  xp  =  *1 

BfCp  =  A21  Cp  =  C-ji  Xf  =  X2 


(1) 


(2) 
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Figure  2-18.  Plant-Filter  Equations 
2.6  DESIGN  EXAMPLE 

These  design  examples  illustrate  the  effectiveness  of  the  filter  accommodation 
algorithm  to  oope  with  residual  modes.  To  keep  the  examples  simple  we  include  one 
mode  from  each  class:  controlled,  suppressed,  and  residual.  In  order  to  guarantee 
residual  mode  excitation  we  make  the  residual  model  have  slopes  identical  to  that  of 
the  controlled  mode.  Thus  we  ensure  that  attenuation  of  the  residual  modal  response 
results  from  the  filter  action  and  not  judicious  actuator  or  sensor  placement.  In 
order  that  the  example  demonstrates  the  compatibility  of  the  filter  accommodation 
technique  and  the  standard  MESS  algorithm,  we  also  include  a  suppressed  mode.  The 
effectiveness  of  the  filter  accommodation  technique  in  promoting  residual  mode 
stability  is  also  judged  by  comparison  with  a  MESS  design  that  does  not  employ 
filter  accommodation.  This  last  step  ensures  that  we  have  not  worked  with  a  benign 


problem.  Eigenvalue  analyses  are  employed  to  judge  the  efficiency  of  the  two  designs. 
The  residual  modes  should  migrate  less  with  the  filter  accommodation.  Time  history 
plots  are  also  used  to  observe  signals  in  the  closed-loop  system.  Less  high  fre¬ 
quency  noise  should  exist  in  the  filter  accommodated  designs.  These  ground  rules 
are  severe  -  especially  the  equal  slope  condition;  however,  the  numerical  results  of 
such  a  design  problem  should  provide  a  clear  indication  of  the  potential  of  the  filter 
accommodation  algorithm. 


(D 


(2) 


(3) 
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Figure  2-19.  Example  Design  Model 

As  shown  in  Figure  2-19,  a  two-mode  design  model  is  employed.  This  model  has 
frequencies  at  one  and  two  radians.  We  choose  to  control  the  one  radian  mode  and 
suppress  control  action  to  the  two  radian  mode.  Thus  the  controlled  state  vector  is 
composed  of  states  x^  and  xg,  and  the  suppressed  state  vector  is  composed  of  states 
x3  and  x4.  The  actuators  and  sensors  are  colocated  therefore  B  =  C  . 
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The  input  matrices  have  the  directions  shown  in  Figure  2-20.  (Bc  =  b2  =  Br  same 
residual  and  controlled  mode  slopes. )  Thus,  BgK  =  ByK  and  the  control  action 
applied  to  the  residual  and  controlled  modes  is  identical.  The  residual  modes  will 
most  certainly  be  excited  and  thus  we  will  be  able  to  check  the  efficiency  of  the 
filters  in  reducing  observation  spillover. 


Figure  2-20,  Modal  Input  Matrix  Row-Vector 
Directions 

One  aspect  of  filter  accommodation  is  the  ability  to  provide  attenuation  which  starts 
deep  within  the  controlled  region  and  still  provide  the  proper  amount  of  control 
authority  to  meet  performance  requirements.  To  demonstrate  this  aspect  of  the 
algorithm  we  choose  a  second-order  filter  that  corners  at  one-tenth  of  a  radian.  The 
controlled  mode  at  one  radian  is  a  decade  away  from  the  filter  corner  frequency  and 
thus  experiences,  when  viewed  from  the  filter,  over  40  decibels  of  attenuation. 

Figure  2-21  illustrates  the  filter  frequency  response. 


FILTER  TRANSFER  FUNCTION 


F(S)  = 


(01X011) 

(S+0.1XS+0.11) 


FILTER  FREQUENCY  RESPONSE 


MODAL  FREQUENCIES 

«c  =  1.0  RAD 
u,  =  2.0  RAD 
wr  =  3.0  RAD 
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Figure  2-21.  Roll-Off  Filter  Characteristics 


2. 6. 1  TIME  HISTORIES.  We  next  apply  a  unit  initial  condition  to  Mode  1  and 
examine  the  transient  response  of  all  signals  around  the  loop.  Figure  2-22  illustrates 
the  response  of  modal  positions:  Mode  1,  the  controlled  mode,  behaves  quite  well  and 
exhibits  the  required  damping  characteristics.  Mode  2,  the  suppressed  mode, 
receives  a  small  amount  of  excitation  (~  10"11)  and  behaves  in  open  loop  fashion. 

Mode  3,  the  residual  mode  is  excited  by  controller  action  (« 10-1)  and  rings  in  open- 
loop  fashion. 


Figure  2-23  illustrates  the  sensor  and  filter  responses.  The  excited  residual  mode 
spills  out  into  the  sensors  and  the  filter  effectively  attenuates  the  residual  mode 
response. 


2. 6. 2  COMPARISON  WITH  A  NONFILTER  ACCOMMODATED  DESIGN.  To  check 
the  efficiency  of  filter  accommodation  to  cope  with  residual  mode  spillover,  a  non¬ 
filter  accommodated  design  employing  the  standard  MESS  algorithm  is  applied  to  the 
design  model.  The  design  parameters  (a-shift  and  suppression  level)  are  identical  for 


RESIDUAL 


0.05 


SENSOR  AND  FILTER  RESPONSES 
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Figure  2-23.  Sensor  and  Filter  Responses 

both  designs.  Eigenvalue  data  for  the  designs  is  presented  in  Tables  2-1  and  2-2. 

Due  to  observation  spillover  the  nonfilter  accommodated  design  is  unstable. 

Instability  occurs  with  the  associated  residual  mode  eigenvalues  of  1.4  ±  j  4.9.  The 
corresponding  residual  mode  eigenvalues  for  the  filter  accommodated  design  are 
-2.2  x  10“2  ±  j  3. 9,  which  indicates  a  small  shift  in  the  real  part  of  the  closed-loop 
root  from  its  open  loop  value  of  -2.0  x  10“2.  Thus,  the  comparison  indicates  that  the 
evaluation  model  was  not  benign  and  demonstrates  the  efficiency  of  the  design 
algorithm  relative  to  this  example.  Figure  2-24  depicts  the  residual  mode  root  migra¬ 
tion  of  the  two  designs  vectorially. 

2. 6. 3  FREQUENCY  RESPONSE  ANALYSIS.  The  filter  accommodation  algorithm 
was  developed  to  provide  bandwidth  control  for  the  LQG  design  process.  With  respect 
to  large  space  system  control,  this  implies  that  all  modes  above  a  certain  frequency 
spectrum  are  gain  stabilized.  To  demonstrate  the  effectiveness  of  the  algorithm,  we 
compare  the  frequency  response  of  the  open-loop  gain  matrix  LQ(s)  for  our  two  design 
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Table  2-1.  Nonfilter-Accommodated  Eigenvalues  (Unstable) 


:igan 

RmI 

Imaginary 

Raaultant 

Anglo 

Damping 

No. 

Part 

Part 

Valua 

(dograoa) 

Ratio 

1 

-13.72 

0. 

13.72 

180.0 

1.000 

2 

1.403 

4.903 

5.100 

74.03 

-.27511  _  , 

Residual 

-.2751 J 

3 

1.403 

-4.903 

5.100 

-74.03 

4 

-1 .9845 

0. 

.9645 

180.0 

1.000 

5 

-.4258 

1.731 

1.783 

103.8 

.2388 

6 

-.4258 

1.731 

1.783 

-103.8 

.2388 

7 

-2.0000E-02 

2.000 

2.000 

90.57 

1  0000E-02]  „ 

„  „  1  Suppressed 

1 .00O0E-02) 

8 

— 2.0000E-02 

-2.000 

2.000 

-90.57 
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Table  2-2.  Filter  Accommodated  Eigenvalues  (Stable) 


Elgon 

Root 

‘~itg ‘~ia~ii 

NnanNnnt 

Angle 

Damping 

No. 

Part 

Part 

Valoa 

(degree*) 

Ratio 

1 

-5.547 

3421 

6  518 

148.3 

.8511 

2 

-5.547 

-3  421 

6.518 

-148.3 

.8511 

3 

-5.962 

0. 

5.962 

180.0 

1.000 

4 

-23242E-02 

3  917 

3.917 

90.34 

5.9333E-03] 

__  _  Residual 

5 

— 2.3242E-02 

-3.917 

3.917 

-90.34 

5.9333E-03J 

6 

-.5654 

1.723 

1.813 

108.2 

.3119 

7 

-.5654 

-1.723 

1.813 

-108.2 

.3119 

8 

-.1000 

0. 

.1000 

180.0 

1.000 

9 

-.1100 

0. 

.1100 

180.0 

1.000 

10 

-.8789 

0. 

.8789 

180.0 

1.000 

11 

-3.227 
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examples.  This  comparison  is  made  for  plants  with  and  without  residual  modes. 
Also  the  frequency  response  of  the  optimal  compensators  are  compared  for  each  de¬ 
sign. 
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Figure  2-24.  Residual  Mode  Root  Migration 

As  both  L0(s)  and  the  compensator  transfer  characteristics  are  characterized  by 
(4  x  4)  matrices,  for  conciseness  of  presentation  we  provide  the  frequency  response 
for  the  one-one  entry  in  each  case.  (The  remaining  matrix  entries  exhibit  similar  be¬ 
havior.  )  To  determine  L0(s)  we  break  the  feedback  loop  at  point  b  of  Figure  2-17, 
excite  the  system  at  actuator  port  number  one  and  measure  the  signal  returning  to 
actuator  one  from  the  compensator.  This  procedure  is  analogous  to  determination  of 
the  return  difference  for  scalar  systems.  In  similar  fashion  we  obtain  the  compensa¬ 
tor  frequency  response  by  breaking  the  feedback  loop  at  both  points  a  and  b  of  Figure 
2-17,  inject  a  signal  at  sensor  port  number  one  and  observe  the  output  signal  that 
would  be  applied  to  actuator  one.  Thus  we  obtain  the  frequency  response  of  the  one- 
one  entry  of  the  compensator  matrix. 
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Figure  2-25  depicts  the  frequency  response  of  the  one-one  entry  of  L0(s)  for  the 
filter  accommodated  design.  The  design  model  of  the  plant  (one  controlled  mode,  one 
suppressed  mode,  and  no  residual  modes)  is  employed  in  the  frequency  response  gen¬ 
eration.  Examination  of  the  plot  reveals  fourth-order  behavior  at  high  frequencies, 
i.e. ,  80  db  roll -off /decade.  This  high  frequency  behavior  can  be  explained  as  follows: 
the  optimal  controller  behaves  as  1/s,  the  optimal  estimator  behaves  as  l/s,  and  the 
filters  behave  as  1/s2  where  s  is  the  complex  Laplace  transform  parameter.  Thus  we 
achieve  fourth-order  high  frequency  behavior  as  our  design  analysis  predicted.  This 
fourth-order  response  is  in  contrast  to  the  standard  optimal  control  algorithm  which 
provides  only  a  second-order  roll-off  at  high  frequencies. 

Figure  2-26  depicts  the  frequency  response  of  the  one-one  entry  of  L0(s)  for  the 
standard  optimal  design.  Again  we  employ  the  design  model  of  the  plant  (no  residual 
modes)  for  the  frequency  response  generation.  As  anticipated,  the  plot  exhibits 
second-order  behavior  at  high  frequencies,  i.e. ,  40  db  roll-off/decade,  which 
corresponds  to  single  integration  behavior  for  both  the  controller  and  the  estimator. 

We  also  note  the  flatness  of  the  response  curve  after  the  resonant  peak  in  Figure  2-26. 
This  response  is  characteristic  of  standard  optimal  control  designs  which  tend  to  be 
high  bandwidth  in  nature.  The  "plateau"  effect  has  been  reduced  in  the  filter  accom¬ 
modated  designs. 

The  resonant  peaks  at  0. 159  Hz  present  in  both  Figures  2-25  and  2-26  correspond  to 
controlled  mode  resonance.  The  suppressed  mode  (0. 318  Hz)  is  not  visible  in  the 
frequency  response  plots.  This  phenomena  occurs  because  the  basic  MESS  optimal 
control  algorithm  employed  in  both  designs  renders  the  suppressed  mode  uncontroll¬ 
able  and  unobservable;  thus  effectively  blocking  excitation  of,  and  transmission  from, 
the  suppressed  mode. 

Additional  response  plots  of  LQ(s)  with  residual  modes  included  are  shown  in  Figure 
2-27  and  Figure  2-28,  for  the  filter  accommodated  design  and  the  standard  optimal  de¬ 
sign,  respectively.  Examination  of  the  residual  mode  resonant  peaks  at  0. 477  Hz 
reveals  the  attenuation  effects  of  the  filter.  The  resonant  peak  corresponding  to  the 
residual  mode  exhibits  a  magnitude  of  3.  5  for  no  filter  accommodation  and  a  magnitude 
of  1. 2  for  the  filter  accommodated  design;  a  difference  of  approximately  8. 6  db.  Thus 
our  filter  accommodated  design  successfully  attenuates  the  residual  mode  response  as 
required. 

We  next  examine  the  frequency  response  of  the  optimal  compensator.  The  compen¬ 
sator  is  defined  as  that  system  which  has  the  sensor  vector  as  its  input  and  the  control 
vector  as  its  output;  i.  e. ,  the  transfer  function  matrix  between  points  "a"  and  "b"  of 


FREQUENCY  (HERTZ) 


Figure  2-25.  Filter  Accommodated  Design.  Loop  Gain  Frequency 

Response  for  L0(s).  (No  Residual  Modes  -  Design  Model) 
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Figure  2-26.  Non-Filter  Accommodated  Design.  Loop  Gain  Frequency 
Response  for  LQ(s).  (No  Residual  Modes  -  Design  Model) 
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Figure  2-27.  Filter  Accommodated  Design.  Loop  Gain  Frequency 

Response  for  L0(s).  (Residual  Mode  Included  -  Evaluation 
Model) 
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Figure  2-17  for  filter  accommodated  design.  For  the  non-filter  accommodated  design, 
the  compensator  is  defined  precisely  as  before  except  the  filters  are  replaced  by 
dynamicless  identity  matrices  and  the  filter  accommodation  matrix  is  a  zero  matrix. 
Figure  2-29  depicts  the  standard  textbook  optimal  compensator. 

Figure  2-30  displays  the  frequency  response  of  the  one-one  entry  for  the  filter  ac¬ 
commodated  compensator.  The  magnitude  of  the  response  initially  remains  flat 
(magnitude  =  0.  6)  and  then  rolls  off  at  60  db/decade;  i.e. ,  third-order  high  frequency 


Figure  2-29.  Standard  Textbook  Optimal  Compensator 

behavior.  This  behavior  is  to  be  expected  as  the  filters  behave  as  1/ s^  at  high  fre¬ 
quencies  and  estimator  behaves  as  1/s.  Thus  the  compensator  behavior  and  the  L0(s) 
behavior  are  in  agreement.  The  additional  matrix  entries  for  the  compensator  re¬ 
sponse,  now  shown  here,  provide  similar  high  frequency  behavior. 
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Figure  2-31  displays  the  frequency  response  of  the  one-one  entry  for  the  standard 
optimal  compensator.  This  response  exhibits  a  peak  at  1. 7  Hz  and  their  rolls  off  at 
20  db/decade;  i.  e. ,  first-order  behavior  at  high  frequencies.  Again  this  high  fre¬ 
quency  behavior  is  expected  for  the  standard  optimal  compensator.  With  respect  to 
the  residual  mode,  the  peak  at  1.7  Hz  represents  phase  lead,  and  is  detrimental  to 
residual  mode  stability;  i.  e. ,  the  system  gain  increases  from  approximately  0. 5  to 
2. 5,  or  14  db.  This  behavior  also  contributes  to  the  flatness  of  the  frequency  re¬ 
sponse  plot  for  Lq(s)  observed  previously  in  Figure  2-26. 

Thus  the  frequency  response  plots  confirm  the  efficiency  of  our  analytical  design 
methodology  and  the  filter  accommodation  algorithm  provides  a-priori  high  frequency 
behavior  as  predicted. 

2.  7  FILTER  ACCOMMODATION  AND  LIAPUNOV  STABILITY 

The  effects  of  the  roll-off  filter  may  be  easily  incorporated  into  a  Liapunov  stability 
analysis.  The  plant  transfer  function  matrix  has  the  form  shown  in  equation  (2)  of 
Figure  2-32.  The  diagonal  nature  of  the  matrix  structure  results  from  the  normal 
mode  representation  of  the  structure.  If  we  perform  the  matrix  multiplication  indi¬ 
cated  in  equation  (2),  we  obtain  the  linear  sums  shown  in  equations  (4)  and  (5).  Thus, 
the  overall  plant  may  be  viewed  as  a  linear  sum  of  the  controlled,  suppressed  and 


residual  plants. 

P(s)=C(sl-A)-1B  (1) 

P(s)=[CcCsCr]  RsI-Ac)1  0  0  I  [Be’ 

0  (Sl-As)’1  0  Bs  (2) 

-  0  0  (sl-ArH  LBr. 

P(s)=Cc(sl-Ac)'1Bc+Cs(sl-As)1Bs+Cr(sl-Ar)'1Br  (3) 

P(s)=Pc(s)+Ps{s)+Pr(s)  (4) 

Attenuate  Pr($)  with  a  filter  such  that  Y2  <<  1 
72Pr(s)=Y2Cr(sl-Ar)‘1Br  (5) 


The  matrix  Cr  has  been  reduced  by  the  attenuation  factor: 

Cr  — ►72Cr; 

15051276*6 

Figure  2-32.  Roll-Off  Filter  Attenuation  and  the  Open-Loop  Plant  Transfer 
Matrix 
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Thus,  if  we  assume  that  for  the  residual  mode  frequency  range  the  filter  has  an 
attenuation  factor  this  factor  may  be  directly  incorporated  into  the  residual 

plant  transfer  function  matrix,  as  shown. 

From  this  perspective,  the  filter  modifies  the  effective  value  of  the  input  and  output 
matrices  by  multiplication  with  a  constant,  i.  e. ,  Cr  —  y  2Cr.  Thus  we  may  incor¬ 
porate  modified  matrix  in  a  Liapunov  analysis  for  stability.  This  analysis  reflects  the 
fact  that  the  filter  attenuates  observation  spillover,  but  does  not  affect  control  spill¬ 
over. 

A  closed-loop  stability  analysis  for  the  closed-loop  system  involves  the  closed-loop 
system  matrix  shown  in  Figure  2-33  where  it  is  assumed  that  we  have  an  estimator 
in  the  loop:  i.  e. ,  xc,  ec,  xs,  and  xr  represent  the  controlled  state  vector,  the  error 
vector  for  the  estimated  controlled  states,  the  suppressed  state  vector,  and  the  re¬ 
sidual  state  vector,  respectively.  A  Liapunov  analysis  can  be  applied  to  this  matrix. 
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Figure  2-33.  Complete  Closed-Loop  System  Matrix  Including 
Residual  Modes  and  Attenuation  Factor 


The  closed-loop  ays  tem  matrix  may  be  diagrammed  as  shown  in  Figure  2-34.  The 
effect  of  the  filter  Is  to  attenuate  the  response  generated  by  the  residual  modes.  For 
zero  attenuation  factor,  the  residual  mode  loop  would  be  completely  open.  As  pre¬ 
viously  noted,  the  residual  modes  are  excited;  control  spillover  is  not  attenuated,  but 
observation  spillover  is  attenuated  by  the  filter  attenuator  factor.  It  is  the  attenuation 
of  observation  spillover  that  counteracts  the  potentially  destabilizing  effect  of  residual 
mode  excitation.  This  unidirectional  decoupling  achieved  by  filter  accommodation  is 
not  as  strong  as  that  achieved  by  the  basic  MESS  algorithm;  however,  the  MESS 
algorithm  relies  on  knowledge  of  the  decoupled  modes.  In  the  face  of  incomplete 
knowledge,  except  perhaps  a  frequency  range,  filter  attenuation  of  the  unknown  modal 
responses  is  the  best  that  can  be  achieved  with  linear  techniques. 


Attenuation 

factor 


Suppressed  Estimator  Controlled 

dynamics  dynamics  dynamics 
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Figure  2-34.  Closed-Loop  System  Dynamics  with  Residual  Modes 

Following  Siljak's  (Reference  22)  decentralized  Liapunov  methodology  we  arrive  at 
the  stability  conditions  indicated  in  Figure  2-35.  The  Liapunov  stability  parameters. 


the  epsilons,  are  direct  functions  of  the  attenuation  factor  gamma.  Thus,  a  con¬ 
servative  requirement  for  the  degree  of  attenuation  required  by  the  roll-off  filter  is 
established. 

The  control  spillover  term  e  21  can  be  expressed  in  standard  matrix  norm  notation, 
as  shown  in  Figure  2-36.  The  Schwarz  inequality  is  then  easily  applied  to  the  matrix 
product  BrK.  This  allows  separation  of  the  norms  for  both  matrices.  Similar  re¬ 
sults  hold  for  observation  spillover.  Thus,  the  bounds  may  be  expressed  as  shown 
by  inequality  (5).  Assumptions  concerning  the  structural  damping  and  modal  slopes 
together  with  gain  levels  allows  calculation  of  the  filter  attenuation  constant,  gamma. 

The  expressions  developed  for  the  stability  analysis  provide  an  analytical  basis  for 
several  heuristic  concepts:  Examination  of  inequality  (5)  reveals: 


(fs^s)  Jr)  t 

- - >  t12t21 

^12*121 


where 

Or 

e12 


€21 


Ql2  Pl2 

*m 


suppressed  damping  (min  value) 
residual  mode  damping  (min  value) 
observer  spillover  term 
^  b2  CJqJ  GCr)  —  7  I  GCr  I 
control  spillover  term 
*2  [2(BrK)T  BrK]  =  V2  II  B^  » 
Liapunov  parameters  (always  positive) 
Maximum  eigenvalue 
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Figure  2-35.  Roll-Off  Filter  Attenuation  and  Liapunov  Stability 
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Figure  2-36.  Spillover  Bounds  as  a  Function  of  System  Parameters 

1)  It  is  always  possible  to  provide  gain  stabilization  for  the  system,  i.  e. ,  as 
gamma  tends  to  zero,  open-loop  conditions  prevail  in  the  LSS  which  is  open- 
loop  stable. 

2)  High  gain  levels  in  both  filter  or  estimator  require  heavy  filtering  (or  a  large 
amount  of  intrinsic  modal  damping)  ensure  stability.  For  a  fixed  amount 

of  modal  damping  and  a  high  gain  system,  the  attenuation  provided  by  the  filter 
must  be  correspondingly  larger  to  ensure  satisfaction  of  the  stability  inequality. 

3)  The  norms  ||  k||  and  ||g||  are  proportional  to  available  actuator  and 
sensor  power.  Thus,  one  can  obtain  stability  bounds  as  functions  of  assumed 
residual  mode  damping  and  residual  mode  slope. 

2. 8  DISTURBANCE  ACCOMMODATION 
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For  those  systems  in  which  modal  damping  of  itself  is  insufficient  to  meet  line  of 
sight  specifications,  and  the  Draper  evaluation  modal  #2  falls  into  this  class  of 
structure,  some  additional  control  technique  must  be  employed  to  meet  performance 
criteria.  Disturbance  accommodation  is  one  technique  that  can  be  employed  to 
counteract  disturbances.  We  choose  disturbance  accommodation  because  it  allows  us 
to  integrate  both  open-loop  and  closed-loop  control  techniques,  i.e. ,  we  can  set  the 
disturbance  accommodation  loop  independently  of  the  damping  augmentation  loop. 

Disturbance  accommodation  falls  basic  into  two  broad  classes:  input  disturbance 
estimation  (E.  J.  Davison,  Reference  23)  and  output  disturbance  estimation  (C.  D. 
Johnson,  References  24  and  25).  The  main  distinction  lies  in  the  site  of  disturbance 
sensing.  One  technique  (input  estimation)  measures  the  disturbance  at  the  disturbance 
site  and  leads  to  open-loop  correction.  The  other  technique  (output  estimation) 
measures  the  effect  of  the  disturbance  at  the  system  output  and  leads  to  closed-loop 
correction.  For  large  distributed  systems,  input  estimation  is  preferred  due  to  the 
long  propagation  time  of  the  system,  i.e. ,  the  disturbance  could  be  operational  for 
long  periods  of  time  before  it  was  sensed  at  the  system  outputs.  In  addition,  the  open- 
loop  correction  will  provide  no  stability  problems  in  the  face  of  large  parameter  var¬ 
iations.  Figure  2-37  illustrates  the  two  concepts.  The  disturbance  canceling  equation 
may  be  easily  derived  as  follows:  First  we  note  that  disturbance  signal  uj  must  cancel 
the  disturbance  effects  for  any  w^.  Thus  we  form  the  equations  shown  in  Figure  2-38. 


THE  DISTURBANCE  STATES  ARE  SENSED  AT  THE  THE  DISTURBANCE  STATES  ARE  SENSED  AT  THE 
DISTURBANCE  INPUT  AND  A  COMPENSATING  SYSTEM  OUTPUT  AND  A  COMPENSATORS  SIGNAL  IS 

SIGNAL  IS  DEVELOPED  THROUGH  K0.  DEVELOPED  THROUGH  THE  DISTURBANCE  ESTIMATOR. 


Figure  2-37.  Disturbance  Accommodation  Techniques 
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ACCOMMODATION  EQUATION:  ud  MUST  CANCEL  wd  FOR  ANY  X6 

84X4  -  BKdXd  =  0  (1) 

OR 

*d  -  BKd  =  0  (2) 
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Figure  2-38.  Input  Disturbance  Accommodation 
Equation  (3)  is  the  disturbance  accommodating  equation. 

As  noted,  input  disturbance  accommodation  is  an  open-loop  technique.  The  accom¬ 
modation  equation  is  static  in  nature.  Another  strength  of  the  procedure  is  that  it 
allows  the  closed-loop  damping  control,  uc,  to  be  implemented  independently  of  the 
disturbance  control. 

The  disturbance  accommodation  loop  is  easily  analyzed  by  transfer  functions.  The 
accommodation  loop  modifies  the  input  matrix  B^.  Thus,  exact  satisfaction  of  the 
accommodation  equation  (3)  will  result  in  zero  disturbance  effects  at  the  system  out¬ 
put,  as  shown  in  Figure  2-39. 

The  exact  solution  of  the  accommodation  equation  requires  inversion  of  the  input 
matrix  B.  In  general,  this  matrix  does  not  have  an  inverse;  however,  when  the 
number  of  actuators  equals  the  number  of  controlled  modes,  an  exact  solution  for  the 
accommodation  equation  may  be  obtained  by  following  the  procedure  outlined  in 
Table  2-3. 
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(1) 


NO  DISTURBANCE  ACCOMMODATION 
V  =  C(SI  -  A)“1B-Xi 

DISTURBANCE  ACCOMMODATION 


Y  =  C(SI-ArMBd-Wd>X|» 

(2) 

THUS 

Bd  -  (8,1  -  BKtf)  =  0 

(3) 

FOR  PARAMETER  VARIATIONS: 

B,)  -  Bd  +  ABd 

W 

B  -  B  +  AB 

(5) 

THUS  ' 

(8^  8Kd)  -  [Bj  +  AB*  -  (B  +  BAlKd  s  [AB^  —  ABK^]  *  0 
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Figure  2-39.  Disturbance  Accommodated  Transfer  Function 

To  clarify  the  procedure  just  outlined,  a  sample  problem  is  presented.  Table  2-4. 

The  example  illustrates  the  form  of  B~*  and  how  exact  satisfaction  of  the  accommoda¬ 
tion  equation  is  obtained. 

_ Table  2-3.  Disturbance  Accommodation  Solution  Procedure _ 

Bd  -  BKo  =  0  (1) 

OR 

*0  =  B-1Bd  (2) 


M  GENERAL,  B  IS  NOT  BIVERT1BIE.  HOWEVER,  WE  WANT  AN  EXACT  SOLUTION  FOR 
Kg  (NO  GENERALIZED  DIVERSE).  THUS,  WE  DO  THE  FOLLOWMG  OPERATIONS: 

1)  FORM  B  BY  0ELETM6  ZERO  ROWS  FROM  B  AND  Bj 

2)  CAUSE  B  TO  BE  SQUARE  BY  MANN8  THE  NUMBER  OF  ACTUATORS  TO  EQUAL  THE 
NUMBER  OF  MOOES  TO  BE  CONTROLLED 

3)  INVERT  B  (WE  PUCE  ACTUATORS  TO  ASSURE  THE  INVERT16RJTY  OF  I) 

4)  INSERT  CORRESPONDING  ZERO  COLUMNS  INTO  jM  FORMMG  B-1 
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As  with  any  open-loop  technique  disturbance  accommodation  may  be  sensitive  to 
parameter  variations.  The  disturbance  accommodation  equation  can  be  desensitized 
against  parameter  variations  by  adding  additional  actuators  such  that  the  original  in¬ 
put  matrix  (which  had  our  special  inverse)  is  imbedded  in  a  larger  matrix.  This 
desensitization  phenomenon  can  be  explained  by  two  effects:  reduced  gain  level  and 
signal  averaging.  If  we  add  additional  actuators,  the  same  general  actuation  level 
may  be  attained  with  lower  gains.  Also,  as  sensor  and  actuator  signals  are  now  re¬ 
ceived  from,  and  applied  to,  a  larger  area  of  the  structure  parameter  variations  tend 
to  "average  out. " 

The  desensitization  effect  can  be  further  enhanced  if  one  has  a-priori  information  of 
areas  of  the  structure  that  are  uncertain:  one  simply  selects  the  weighting  matrix  R 
to  cause  small  gains  for  the  troublesome  areas.  (This  weighting-matrix  procedure 
was  irot  used  for  the  Draper  evaluation  structure  because  it  violated  the  ground  rules. ) 

To  derive  the  desensitized  disturbance  gains  we  assume  the  number  of  actuators 
exceeds  the  number  of  controlled  modes,  but  in  placing  our  actuators  we  retain  a 
square  invertible  matrix,  Bj  (n  x  n)  which  has  rank  equal  to  the  number  of  controlled 
modes.  Thus  is  imbedded  in  the  larger  actuator  matrix,  B  (2n  x  r),  and  we  now 
proceed  to  solve  this  augmented  filter  accommodation  equation  for  the  disturbance 
gain  matrix  Kd.  The  procedure  is  detailed  in  Table  2-5;  where  we  employ  a  specific 
optimal  control  format  thereby  generalizing  a  procedure  given  by  Kwaakemak  and 
Silvan,  p.  262,  Reference  17;  i.e. ,  Pj2  is  a  submatrix  of  the  Riccati  equation.  We 
also  note  that  the  solution  obtained  for  Kd  is  exact. 

2.9  INPUT-OUTPUT  DECENTRALISATION 

We  now  take  advantage  of  the  many  actuators  required  for  disturbance  accommodation 
to  design  a  completely  decoupled  control  law  such  that  each  mode  is  controlled  by  a 
separate  gain  elemeit  (Reference  22).  As  both  B  and  C  are  invertable  in  the  sense  in¬ 
dicated,  and  we  have  an  actuator/sensor  pair  for  each  mode,  no  estimator  is  required. 
This  leads  to  a  simple  closed-loop  controller  of  increased  robustness. 

Figure  2-4  illustrates  the  concept  where  we  show  the  inversion  of  the  input  and  output 
matrix  in  conjunction  with  a  diagonal  feedback  gain  matrix;  i.  e. ,  the  actual  system 
gains  are  given  by  the  matrix  (-B"1  K  C-1). 

The  effects  of  input-out  decentralization  are  easily  analyzed  in  state  variable  form. 
The  matrix  inversions  for  B  and  C  effectively  eliminate  these  matrices  from  the 
dynamic  equations.  If  the  gain  matrix,  K,  is  diagonal  then  each  mode  is  separately 
controlled,  i.e. ,  we  have  completely  decoupled  control.  Figure  2-41  illustrates  the 
procedure. 
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Table  2-4.  Sample  Calculation  of  the  Disturbance 


i 
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Table  2-5.  Desensitizing  Effect  of  Additional  Actuators 

LET  THE  NUMBER  OF  ACTUATORS  EXCEED  THE  NUMBER  OF  MOOES,  BUT  WE  RET  AM 
OUR  SQUARE  MVERTWLE  MATRIX 


B-|(N  x  N)  -  B(2N  x  R)  =  [B-)  B2];  B— 1  EXISTS,  R  >N  (1) 

NOW  SOLVE  THE  ACCOMMODATION  EQUATION  BY  GENERALIZED  INVERSE  AND 
SPECffC  OPTIMAL  CONTROL  FORMAT 


OR 

Kj  =  R-1BTP12 

(2) 

AND 

BK0  =  BR“lBTP12  =  B,j 

(3) 

THUS 

Pl2=  |BR-“1BTr1Brf 

(4) 

LET  (R  =  0 

Kd=  R-1BT[BR“lBT)“1Bd 

(5) 

AND 

Ko= 

(6) 

B,,  =  8X4  (EXACT  SOLUTION) 

(7) 

264.7604 


2. 10  BANDWIDTH  REDUCTION  AND  DISTURBANCE  ACCOMMODATION 

The  disturbance  control,  ly,  and  the  closed-loop  control,  uc,  are  independent  of 
each  other.  Thus,  u£  can  be  designed  by  any  means  available  to  the  designer, 
optimal,  suboptimal  or  classical  means.  Of  special  interest  is  the  design  of  low 
bandwidth  systems  via  input/output  decentralization.  Such  decoupled  systems  will 
provide  greater  robustness  in  the  face  of  severe  parameter  variations. 

Transfer  function  analysis  also  reflects  this  independence  of  uc  and  u^.  The  disturb¬ 
ance  control  modifies  the  input  matrix,  B^.  The  closed-loop  control  modifies  the 
plant  matrix,  A.  Both  modifications  appear  as  independent  multiplicative  factors  in 
the  transfer  function,  as  shown  in  Figure  2-42. 
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DECENTRALIZED  INPUT 


'DECENTRALIZED  OUTPUT 


264  760-7 

Figure  2-40.  Input-Output  Decentralization  for  Decoupled  Control 


X  =  Ax  +  Bu  (1) 

7=  Cx  (2) 

Uc  =  -B-^C-’T  =  -B-1KC-1CX  =  -B_1KX  (3) 

THUS 

X  =  Ax  +  BI-B-’KX)  (4) 

X  =  Ax  -  KX  (51 


ASSUMING  A  DIAGONAL  GAM  MATRIX  WELDS  DECOUPLED  CONTROL 
(EACH  GAIN  AFFECTS  ONLY  ONE  MOOE) 
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Figure  2-41.  Effect  of  Input-Output  Decentralization 


CLOSED-LOOP  TRANSFER  FUNCTION 


FREQUENCY 

RESPONSE 


Y  =  C(SI  -  Acl)-1  (Bu  -  BKd)Xd 

_ )  t _ 


DISTURBANCE 

CONTROL 
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Figure  2-42.  Simultaneous  Bandyridth  Reduction  and  Disturbance 
Accommodation 
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SECTION  3 

ANALYSIS  AND  SIMULATION 


The  purpose  of  the  analysis  and  simulation  task  was  to  evaluate  the  applicability  of 
the  theory  to  complex  structure.  Since  erroneous  conclusions  can  be  drawn  from 
evaluations  using  too -simple  models,  the  complexity  of  the  structural  model  is  an 
essential  feature  of  this  task. 

3. 1  DRAPER  MODEL  #2 

In  order  to  provide  a  realistic  testbed  for  control  system  evaluation,  the  Charles 
Stark  Draper  Laboratory  developed  the  system  model  shown  in  Figure  3-1  and 
described  in  detail  in  Reference  3.  This  model  consists  of  a  flexible  optical  support 
structure  and  an  isolated  equipment  section  which  contains  solar  panels.  Optical 
support  structure  consists  of  the  upper  mirror  support  truss,  lower  mirror  support 
truss,  and  the  metering  truss  which  maintains  mirror  separation.  The  finite  element 
model  is  shown  in  Figure  3-2.  The  structure  is  configured  as  a  truss,  but  it  is 
assumed  that  all  joints  allow  a  full  moment  connection.  Thus,  both  bending  and  axial 
stiffness  are  included  in  the  model  for  all  structural  members.  To  evaluate  system 
performance  two  classes  of  disturbances  (parameter  variations  and  oscillatory  equip¬ 
ment  vibrations)  are  included  for  application  to  the  model. 

3.1.1  PARAMETER  VARIATIONS.  To  assess  the  sensitivity  of  the  control  system 
to  modeling  errors  and  manufacturing  tolerances  two  sets  of  physical  parameter 
variations  are  employed.  These  variations  are  summarized  in  Table  3-1.  The 
changes  are  confined  to  the  optical  support  structure  and  involve  mass  distribution 
changes  on  the  upper  optical  support  truss  and  stiffness  changes  in  the  metering  truss. 
These  parameters  are  chosen  based  on  their  impact  on  LOS  error  in  several  critical 
modes.  The  amount  of  variation  is  driven  by  the  desire  to  vary  the  natural  frequen¬ 
cies  by  10%. 

Dynamic  analyses  were  conducted  on  both  perturbed  models  using  NASTRAN.  The 
resulting  natural  frequencies  are  given  in  Table  3-2  along  with  those  of  the  original 
model.  The  results  show  that  both  sets  of  parameter  variations  had  more  impact  on 
the  higher  frequency  modes  because  the  upper  support  truss  and  the  metering  truss 
contain  a  higher  percentage  of  the  total  energy  in  these  modes  than  in  the  lower  modes. 

3.1.2  SYSTEM  DISTURBANCES  AND  OPTICAL  TOLERANCES.  The  disturbance 
model  consists  of  two  sinusoidal  forces  which  are  applied  to  a  point  on  the  optical 
support  structure  and  a  point  on  the  equipment  section.  These  forces  simulate  on¬ 
board  vibrating  equipment.  It  is  assumed  that  both  forces  will  act  simultaneously. 


Figure  3- 
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The  larger  of  the  two  is  located  at  node  46  on  the  equipment  section.  A  smaller 
amplitude  disturbance  is  located  at  node  37  on  the  optical  support  truss.  Both  forces 
act  in  the  z-direction  and  are  assumed  to  act  independently.  The  amplitude,  fre¬ 
quency,  and  direction  of  each  are  given  by  Figure  3-3. 

Damping  is  assumed  to  be  0. 1%  in  all  modes  except  the  isolator  modes  (7,  8,  11,  12, 
13,  16)  which  are  assumed  to  have  50.%  damping  due  to  the  isolator.  The  tolerances 
for  system  line  of  sight  error  and  defocus  are 

-•6 

LOS  =  1. 0  x  10  radians 
-4 

Defocus  =  5.0  x  10  meters 

3. 1. 3  PROBLEM  DESCRIPTION.  The  control  problem  as  presented  by  Draper 
Model  Number  2  may  now  be  stated  as:  Design  a  control  system  for  the  nominal 
structure  which  counteracts  the  two  sinusoidal  disturbances  and  meets  the  LOS 
requirements.  Evaluate  system  performance  in  the  face  of  the  enumerated  parameter 
variations. 

3. 2  INITIAL  PROBLEM  EVALUATION 

Rather  than  proceeding  directly  to  active  system  design,  an  initial  evaluation  of  the 
problem  was  made  as  follows: 

1)  Conduct  an  LOS  analysis  and  rank  modes  according  to  LOS  contribution. 

2)  Arbitrarily  add  20%  damping  to  all  modes  for  test  runs  to  evaluate  potential  role 
of  damping. 

3)  Use  actuator  placement  algorithm  to  provide  input  matrices  of  the  proper  rank. 

4)  Initially  accommodate  all  modes  and  then  delete  mode  by  mode  according  to  LOS 
ranking. 

5)  Add  actuators  to  provide  desensitizatton. 

6)  Analyze  perturbed  systems  to  determine  the  additional  amount  of  attenuation 
required  to  meet  specifications. 

7)  Design  reduced  bandwidth  system. 
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This  evaluation  resulted  in  the  following  conclusions: 

1)  Modal  damping  will  not  provide  disturbance  control  for  this  system. 

2)  Disturbance  accommodation  requires  feedforward  control  of  43  modes.  (Mode 
50  is  residual.) 

3)  Disturbance  accommodation  and  system  bandwidth  reduction  will  be  required  to 
meet  specs  with  parameter  variation,  but  bandwidth  reduction  is  so  large  that  a 
high  gain  system  results  with  the  usual  sensitivity  problems. 

The  negligible  effect  of  damping  on  the  performance  characteristics  may  be 
explained  as  follows:  The  disturbance  frequencies  at  5  and  10  hertz  lie  well  below 
the  resonant  frequencies  of  many  modes.  Thus  for  this  situation  the  modal  transfer 
functions  respond  as  essentially  constant  gains  and  the  disturbances  are  transmitted 
to  the  line  of  sight  scaled  by  the  low  frequency  modal  gain,  i.e. ,  the  modes  act  as 
rigid  links  to  the  LOS.  This  situation  is  illustrated  in  Figure  3-4.  One  obvious  way 
to  correct  this  situation  is  the  lowering  of  the  low  frequency  modal  gain  by  closed- 
loop  means.  This  may  be  accomplished  by  acceleration  feedback  integrated  into  our 
decoupled  control  scheme. 

3.3  DISTURBANCE  ACCOMMODATED  DESIGNS 

Four  different  designs  for  the  Draper  evaluation  structure  are  presented.  Design  I 
involves  standard  disturbance  accommodation  and  employs  43  actuators.  Design  II 
involves  standard  disturbance  accommodation  and  employs  153  actuators.  The 
additional  actuators  are  added  to  provide  desensitization  for  the  disturbance  accom¬ 
modation  algorithm.  Damping  is  injected  in  open-loop  fashion  so  as  to  check  the 
efficiency  of  disturbance  accommodation  without  adding  the  additional  variable  of 
closed-loop  behavior.  Design  HI  incorporates  closed-loop  bandwidth  reduction  in 
addition  to  basic  disturbance  accommodation.  Finally,  Design  I-A  (an  iteration  on 
the  first  design)  provides  20%  damping  via  closed-loop  decoupled  control.  The  final 
design  remains  stable,  but  is  out  of  spec  for  the  perturbed  models.  Table  3-3 
summarizes  the  different  designs  and  results. 


Numerical  evaluation  data  for  each  of  the  designs  is  shown  in  Table  3-4.  Examina¬ 
tion  of  the  table  reveals  that  each  of  the  designs  easily  meets  specifications  by 
several  orders  of  magnitude:  i.  e. ,  nominal  LOS  values  for  Designs  I,  IA,  and  EH 
are  on  the  order  of  10"14  radians.  The  LOS  values  for  Design  EH,  the  low  frequency 
design,  are  on  the  order  of  10~18  radians.  However,  all  designs  are  out  of  specifi¬ 
cation  for  the  parameter  variations.  The  best  performance,  with  respect  to 
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Table  3-3.  Disturbance  Accommodation  Designs 


VELOCITY  FEEDBACK  20%  DAMPING  PERTURBED  MODEL  II:  OUT  OF  SPECIFICAT 

43  ACTUATOR/SENSOR  PAIRS  PERTURBED  MODEL  IV:  OUT  OF  SPECIRCA 

-  43  CONTROLLED  MOOES  ALL  MODELS  STABLE 

-  1  RESIDUAL  MODE  (MODE  SO) 


Table  3-4.  Numerical  Evaluation  of  Design  Data 


r  .  :  i 


parameter  variation,  is  provided  by  Design  n.  The  LOS  error  is  on  the  order  of 
10“6  radian. 

3.4  SYSTEM  DESIGN  COMMENTS 

As  a  result  of  the  analysis  and  simulation  activities  using  the  Draper  optical  system 
model,  it  became  apparent  that  the  system  design  approach  of  designing  the  structure 
first  and  then  adding  the  controls  as  an  "after  thought"  has  severe  limitations  for 
precision  space  systems.  The  best  approach  would  be  one  where  the  structure  and 
the  control  system  are  designed  concurrently  so  to  give  an  optimum  combination  of 
structure  and  control.  However,  since  the  technology  for  the  best  approach  Is  not 
presently  available,  it  would  appear  appropriate  to  at  least  re-examine  Ihe  structure 
for  possible  improvements  where  it  is  found  that  modal  damping  will  not  provide  a 
sufficiently  quiet  optical  environment  to  meet  requirements. 


Ihe  Draper  evaluation  structure  exemplifies  this  viewpoint.  The  disturbances  at 
nodes  37  and  46  are  strongly  coupled  to  the  line  of  sight.  Several  simple  and  econom¬ 
ical  structural  design  changes  can  attenuate  the  disturbance  transmission  character¬ 
istics  such  that  an  active  control  system  can  then  be  appropriately  employed  to  meet 
mission  requirements.  This  conclusion  was  reached  after  observing  the  effect  of  the 
disturbance  accommodation  algorithm:  i.e. ,  it  appears  to  be  attempting  to  redesign 
the  structure,  via  input  matrix  manipulation,  as  viewed  from  the  disturbance  ports. 
When  this  situation  arises,  it  appears  intuitively  more  efficient  to  attack  the  problem 
at  its  primary  source,  the  structural  design  process,  rather  than  indirectly  through 
the  control  system  in  ad  hoc  fashion. 

The  key  phrase  in  the  above  discussion  is  I 'modal  damping".  Again,  it  becomes 
evident  during  the  controller  design  process  that  control  techniques  designed  primar¬ 
ily  for  modal  damping  do  not  provide  an  adequate  design  basis  for  those  classes  of 
systems  that  do  not  respond  to  damping  techniques.  Thus,  the  Draper  model  calls 
attention  to  the  fact  that  multivariable  control  algorithms,  which  incorporate  dis¬ 
turbance  accommodation  to  complement  modal  damping,  may  be  required  for  certain 
classes  of  LSS  control  problems.  This  poses  an  interesting  area  for  future 
research  especially  in  the  arena  of  LQG  design. 
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SECTION  4 


DEMONSTRATION  PLANNING 


Analyses  and  computer  simulations,  while  extremely  valuable  in  Investigating  the 
applicability  of  control  theory,  do  not  provide  the  same  confidence  as  a  comprehen¬ 
sive  demonstration  and  evaluation  test  using  real  hardware.  However,  the  question 
of  the  comprehensiveness  of  the  demonstration  must  be  considered.  The  particular 
structural  configuration  used  in  the  test  could  represent  an  easy  problem  for  the 
theory  and  false  confidence  in  the  ability  of  the  theory  to  handle  other  configurations 
might  result. 

The  purpose  of  the  Demonstration  Planning  task  was  to  identify  the  various  structural 
elements  which  should  be  included  in  a  comprehensive  test  article  and  then  design 
such  a  test  article.  Control  actuator  installation  and  structural  test  article  suspen¬ 
sion  are  also  included  in  the  task.  The  following  were  established  as  goals: 

a.  Configuration  shall  be  representative  of  typical,  projected  space  structure 
systems. 

b.  Elements  of  the  test  structure  shall  be  representative  of  projected  space 
structure  elements. 

c.  Test  structure  shall  exhibit  full  spectrum  of  dynamic  response  and  its  dynamic 
characteristics  shall  be  sufficiently  severe  and  complex  that  the  full  capability 
of  the  control  system  shall  be  adequately  demonstrated. 

d.  Test  structure  shall  have  provisions  for  installation  of  devices  for  propagating, 
sensing  and  controlling  dynamic  disturbance. 

4. 1  SURVEY  OF  SPACE  STRUCTURE  CONFIGURATIONS  AND  ELEMENTS 

A  survey  of  typical,  projected  space  structure  systems  reveals  that  structural 
elements  fall  into  four  general  categories: 

a.  "Plates"  -  for  providing  uniform,  stable  support  to  functional  surfaces  such  as 
mirrors  and  R.F.  reflectors.  For  large  space  systems  these  "plates"  take 
the  form  of  planar  trusses. 

b.  Compression  Members  -  for  reacting  axial  and/or  bending  loads,  i.e. ,  struts 
and  beams,  (pin  ended  or  built-in,  and  cantilevering). 
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c.  Tension  Members  -  for  reacting  axial  tension  loads  only,  l.  e. ,  cords  and 
cables. 

d.  Biaxial  Tension  Members  -  provide  flat,  lightweight,  stable  surfaces;  i.e. , 
membranes. 

Concepts  surveyed  are  listed  in  Table  4-1  and  briefly  described  below. 

a.  HAL/O/WALRUS  (Struts  and  Ties) 

A  structural  system  with  three  major  planar  trusses  stably  located  by  a  fully 
triangulated  system  of  axially  loaded  struts  and  ties. 

b.  HALO/WALRUS  (Cantilevering  beams) 

An  alternative  structural  system  consisting  of  cantilevering  fixed-end  beams 
acting  as  strut  columns;  stabilized  in  torsion  by  triangulated  tension  ties. 

c.  HALO/WALRUS  (Monocoque  shells) 

Another  structural  approach  using  planar  truss  structures  wrapped  through  360° 
to  form  closed  section  "shells"  which  react  loads  in  bending  and  in  shear. 

d.  GEO-STATIONARY  PLATFORM  (Cantilevering  beams) 

Consists  of  six  similar  truss  beams  cantilevering  from  a  common  hub  structure 
with  payload  systems  at  their  tips. 

e.  "ADOPT"  CASSEGRAIN  OPTICAL  CONCENTRATOR  (Trusses  and  beams) 

The  primary  optical  surface  is  supported  by  a  three-tier  system  of  planar 
trusses  and  cantilever  beams.  The  secondary  optical  surface  is  supported  by 
converging  beams  cantilevered  from  the  periphery  of  the  main  "primary"  truss. 
A  large,  high  frequency  planar  truss  provides  "foundation"  structure  for  support 
of  the  primary  optical  surface.  The  interface  between  the  planar  truss  and  the 
primary  optical  surface  Involves  second  and  third  tiers  of  structure.  The 
second  tier  is  an  array  of  truss  "petals, "  and  the  third  is  a  series  of  small 
planar  trusses  that  directly  support  the  optical  panels.  The  "foundation"  planar 
truss  also  provides  support  for  the  beams  that  support  the  system's  secondary 
optical  element. 

f.  FLEXIBLE  FLYER  (Trusses,  beams,  struts  &  ties) 

This  structural  concept  consists  of  a  primary  planar  truss  with  two  systems  of 
triangulated  struts  and  ties  mounted  on  its  periphery,  one  to  support  the  second¬ 
ary  mirror  and  the  other  the  equipment  module.  "Flying"  struts  are  added  to 
reduce  the  effective  length  of  the  major  strut/beams.  All  struts  are  of  the  open 
truss  type. 


g.  MILLIMETERWAVE  RADIOMETER  -  M.R.S.  (Struts,  ties  and  membranes) 

This  750-meter  diameter  system  consists  of  a  torus  supporting  a  "geodesic 
dome"  truss,  which  in  turn,  supports  electrostatically  shaped  reflector 
membranes.  Two  pin-ended  legs  support  the  multiple  feed  assembly.  Each  leg 
consists  of  a  multiplicity  of  pin-ended  struts  and  cross  braced  tension  ties. 

Four  larger  tension  ties,  installed  as  guy  lines,  complete  the  feed  assembly 
support  system. 

h.  "MAYPOLE"  HOOP  AND  COLUMN  (Struts,  ties  and  membrane) 

An  efficient  structural  arrangement  where  diametric  tension  ties  stabilize  the 
compression  hoop.  For  axial  stability  a  central  hub  separates  the  ties  fore  and 
aft  into  conical  orientation.  The  hub  extends,  telescopically,  both  in  the  fore 
and  aft  directions.  The  forward  extension  acts  as  a  cantilever  beam  supporting 
the  R.F.  feed  at  its  tip.  A  membrane,  stretched  within  the  hoop,  is  pulled  into 
near  paraboloidal  shape  by  numerous  tie  lines  anchored  to  the  tip  of  the  aft  hub 
extension. 

i.  LARGE,  LOW-FREQUENCY  TELESCOPE-LOFT  (Struts,  ties  and  membrane) 

A  system  of  tie  lines  and  membranes  held  in  stable  shape  by  centrifugal  force. 
The  axial  components  of  the  centrifugal  forces  are  rea  cted  by  a  single,  com¬ 
pression  element  (mast)  that  is  coincidental  with  the  axis  of  rotation. 

j.  MICROWAVE  POWER  TRANSMISSION  SYSTEM  -  MPTS  (Truss,  struts  and 
membrane) 

This  structure  is  a  NASA  Johnson  Space  Center  design  for  a  Solar  Power 
Satellite  transmitting  antenna.  It  is  a  planar  truss  a  kilometer  in  diameter  and 
consists  of  correspondingly  large,  component,  pin-ended  struts.  Such  struts 
may  themselves  be  of  open  truss  construction. 

Conclusions  to  be  drawn  from  the  survey  are  that  planar  trusses  are  required  in 
almost  all  concepts,  followed  by  pin-ended  struts,  tension  members  (ties),  biaxial 
tension  members  (membranes),  cantilever  beams  and  fixed  ended  strut/beams,  in 
that  descending  order  of  utilization. 

4. 2  TEST  ARTICLE  DESCRIPTION 

Having  identified  the  required  structural  elements,  the  next  step  is  to  incorporate 
them  into  a  structural  test  article.  A  strong  candidate  for  representing  the  planar 
truss  category  is  the  Tetrahedral  Truss  concept.  The  tetrahedral  truss  geometry 
is  recognized  as  a  highly  efficient  space  structure  configuration,  has  been  thoroughly 
characterized,  and  is  likely  to  receive  wide  application.  The  structure  of  the  MPTS 
concept  is  a  very  large  example  of  this  configuration. 
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Table  4-2  presents  dynamic  characteristics  of  a  much  smaller,  eight  bay  tetrahedral 
truss.  The  thirty  degree  orientation  of  its  intermediate  (core)  struts  results  in 
shallower  structural  depth  and  lower  stiffness  than  the  45°  or  60°  orientations. 
Although  the  frequencies  may  change  with  size,  depth,  and  material,  the  modal 
density  should  still  be  typical  for  a  planar  truss.  Thus,  the  close  grouping  of  certain 
modes  (e.  g. ,  7-8,  11-13,  19-20,  21-23,  31-32,  35-36,  38-39)  shows  that  the  tetra¬ 
hedral  truss  has  the  potential  for  providing  a  satisfactorily  complex  control  problem. 

A  tetrahedral  truss  design  is  included  in  candidate  demonstration  hardware  designs 
presented  below.  As  shown  in  Figure  4-1  it  is  an  eight  bay  version,  3. 0  meters 
(10  feet)  in  size,  and  is  candidate  for  representing  the  planar  truss  (plate)  category 
of  space  structure  element.  Material,  sizes  and  proportions  are  selected  to  result 
in  a  first  mode  frequency  below  50  Hz  when  22. 7  Kg  (50  Lb)  of  parasitic  mass  are 
equally  distributed  over  the  face  of  the  structure. 

Due  to  its  high  degree  of  structural  redundancy,  individual  strut  elements  of  the  truss 
may  be  removed  without  loss  of  structural  integrity.  Such  removal,  however,  does 
reduce  the  stiffness  of  die  structure.  It  is  therefore  feasible  to  adjust  and  modify  the 
dynamic  characteristics  of  the  truss  by  discrete  removal  of  appropriate  members 
from  any  location  in  the  structure.  Taken  to  the  extreme,  this  capability  enables 
extensive  modification  to  the  structure,  as  shown  in  Figure  4-2. 

Typical  structural  detail  is  shown  in  Figure  4-3.  The  component  elements  (tubes, 
sleeves  and  discs)  are  common  throughout.  The  structure  consists  of  thin-walled 
brass  tubing,  sleeves  (serving  as  tube  splices)  and  discs  (serving  as  joint  splices  at 
the  structural  node  points).  Ease  of  fabrication  and  of  subsequent  modification  is 
the  driver  in  selection  of  brass.  Silver  soldering  is  selected  for  the  structural 
joints.  Soft  silver  solder  should  be  used  for  the  sleeved  joints  connecting  the  tubular 
elements  to  the  "spiders, "  to  facilitate  possible  subsequent  removal.  The  use  of 
soldered  joints  throughout  minimizes  the  possibility  of  friction  damping  (slip)  at  the 
joints. 

It  is  considered  that  any  alternative  joining  technique  based  on  mechanical  clamping, 
bolting  or  pinning  would  be  cumbersome,  unreliable,  and  expensive  and  may  result 
in  erratic  damping.  There  are  840  sfeeved  joints  and  the  time  and  cost  aspects  of 
selecting,  acquiring,  installing  and  qualifying  that  number  of  mechanical  clamping 
devices  is  likely  to  be  unacceptable. 

Tension  ties  and  cantilever  beams  are  represented  as  shown  in  Figure  4-4.  The 
three  dimensional  truss  beam  is  rigidly  mounted  on  the  planar  truss  at  three  inter¬ 
face  node  points.  Two  guy  lines  (tension  ties)  connect  the  tip  of  the  cantilever  to  • 
points  on  the  periphery  of  the  '  .lanar  truss.  By  varying  the  pre-tension  in  the  ties. 
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Table  4-2.  Typical  Planar  Truss  Modes 
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Figure  4-1.  Test  Structure  Configuration 


Figure  4-4.  Test  Article  with  Beam  and  Tension  Ties 


and  by  removing  the  ties,  a  range  of  overall  dynamic  performance  is  obtainable.  An 
alternative,  "braced  mast"  arrangement  shown  in  Figure  4-5  represents  an  assembly 
of  planar  truss,  pin-ended  strut,  and  three  guy  lines.  Modification  of  dynamic  per¬ 
formance  is  readily  accomplished  by  varying  length  of  the  strut  and/or  tension  in  the 
guy  lines. 

Figure  4-6  illustrates  a  possible  structural  set-up  for  demonstration  testing  that 
utilizes  the  configuration  shown  in  Figure  4-5.  The  structural  model  is  spring 
mounted  on  a  dense  block  representing  the  parent  spacecraft  mass.  Integral  with  the 
block  are  two  tubular  rods,  each  five  feet  long,  representing  low  stiffness  solar 
panel  assemblies.  Cylindrical  weights  are  mounted  on  the  rods  to  simulate  the  mass 
of  the  solar  cells.  The  location  of  these  weights  on  the  rods  is  variable  to  enable 
adjustment  of  the  dynamic  response.  The  entire  assembly  is  suspended  by  means  of 
one  or  two  tension  ties  attached  to  a  point  (or  points)  just  above  the  c.g.  of  the 
assembly.  The  single,  bifurcated  suspension  shown  provides  a  degree  of  pendulum 
stability  in  pitch.  The  degree  of  pendulum  stability  in  the  roll  direction  is  dependent 
upon  the  height  of  the  point  of  bifurcation  above  the  c.g.  Stability  in  yaw  is  dependent 
on  the  torsion  pendulum  stiffness  of  the  suspension  tie.  Yaw  stability  can  be  much 
increased  by  the  use  of  two  parallel  suspension  ties,  one  from  each  attachment  point. 

4.  3  ACTUATOR  INSTALLATION 

Torque  wheels  were  chosen  for  the  control  actuator  installation  studies  since  they 
are  the  only  known  device  which  can  be  used  at  both  high  and  low  frequencies.  Con¬ 
trol  torque  is  generated  by  accelerating  the  rotor  inertia  and  reacting  die  resulting 
equal  but  opposite  torque  back  into  the  structure. 

One  possible  torque  wheel  installation  is  shown  in  Figure  4-7.  It  can  be  installed 
on  any  group  of  three  nodes  anywhere  on  the  back  or  front  face  of  the  structure.  This 
"triple  node"  mount  design  permits  application  of  bending  loadings  to  the  structure 
without  introduction  of  bending  loads  in  the  individual  tubular  elements  of  the  structure. 
The  complete  installation  bolts  directly  to  three  node  points  of  the  demonstration 
structure.  Torque  motor  loads  are  reacted  by  normal  loads  at  the  three  node  points 
and  the  tubular  elements  experience  only  axial  loading.  The  support  provisions 
consist  of  three  bifurcated  legs  providing,  in  effect,  tripod  support  of  both  the  upper 
and  lower  attachment  points  of  the  actuator  assembly.  The  trunnion-type  mounting 
of  the  actuator  assembly  permits  it  to  be  positioned  in  any  desired  orientation  about 
the  trunnion  axis.  The  ring  plate,  supporting  the  stator,  is  provided  with  two  attach¬ 
ment  lugs  180  degrees  apart. 

An  alternative  torque  wheel  installation  is  shown  in  Figure  4-8  involving  only  one  node 
of  the  structure.  Torque  motor  loads  are  reacted  by  bending  loads  in  the  nine  tubular 
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•gure  4-6.  Test  Installation 


Single  Node  Actuator  Installation 


structural  elements  that  converge  at  (he  node  point.  An  adapter  housing,  consisting 
of  an  annular  plate,  a  saddle  bracket,  and  a  circular  base  provides  a  rigid  interface 
between  the  stator  (outer  component  of  the  torque  motor)  and  the  structure  at  the 
node  point.  The  circular  base  plate  fits  over  the  node  assembly,  and  its  circular 
pattern  of  twelve  holes  enables  clamped  attachment  to  the  tubular  elements  of  the 
structure  by  means  of  six  ”U"  bolts.  Since  the  twelve  holes  are  equally  spaced,  the 
torque  motor  assembly  may  be  clocked  in  30°  increments  to  any  preferred  orienta¬ 
tion. 
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